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ABSTRACT 


An  analytical  procedure  is  developed  for  jet  flapped  wings,  in 
which  three-dimensional  features  of  the  trailing  vortex  system  are  repre¬ 
sented.  More  conservative  induced  drag  is  predicted  than  by  analyses  based 
on  the  traditional  planar  vortex  system.  Illustrative  examples  are  given 
showing  that  this  method  gives  better  agreement  with  measured  drag  than  the 
planar -vortex  methods. 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 


1.1  Introduction 

The  U.S.  Air  Force’s  need  for  modernization  of  its  Tactical 
Airlift  capability  led  to  establishment  of  the  Tactical  Airlift 
Technology  Advanced  Development  Program  (TAT-ADP) ,  contributing  to 
the  technology  base  for  development  of  an  Advanced  Medium  STOL 
Transport  (AMST) . 

The  AMST  must  be  capable  of  handling  substantial  payloads 
and  using  airfields  considerably  shorter  than  those  required  by  large 
tactical  transports  now  in  the  Air  Force  inventory.  If  this  short 
field  requirement  is  to  be  met  without  unduly  compromising  aircraft 
speed,  economy,  and  ride  quality,  an  advanced- techno logy  powered-lift 
concept  will  be  required. 

The  STOL  Tactical  Aircraft  Investigation  (STAI)  is  a  major 
part  of  the  TAT-ADP,  and  comprises  studies  of  the  aerodynamics  and 
flight  control  technology  of  powered-lift  systems  under  consideration 
for  use  on  the  AMST.  Under  the  STOL-TAI,  The  Boeing  Company  was  awarded 
Contract  No.  F33615-71-C-1757  by  the  USAF  Flight  Dynamics  Laboratory  to 
conduct  investigations  of  the  technology  of  the  vectored-thrust  and 
Internally  blown  jet  flap  powered-lift  concepts.  These  investigations 
included : 

o  Aerodynamic  analysis  and  wind  tunnel  testing 
o  Configuration  studies 

o  Control  system  design,  analysis,  and  simulation 

This  report  presents  the  results  of  an  analytical  investiga¬ 
tion  of  the  aerodynamics  of  jet  flapped  wings,  with  emphasis  on  drag 
behavior  in  conditions  where  the  classical  assumption  of  an  essentially 
planar  system  of  wing,  jet,  and  trailing  vortex  system  no  longer  gives 
satisfactory  results.  The  analysis  method  was  intended  to  apply  not 
only  to  wings  with  "internally  blown"  jet  flaps,  but  also  to  wings 
with  externally  blown  jet  flaps  or  with  upper  surface  blowing.  J 

1.2  Summary 

An  analytical  procedure  was  developed  to  calculate  the  span- 
wise  load  distribution,  forces,  and  moments  on  wings  embodying  the  jet 
flap  powered,  lift  concept.  The  procedure  is  programmed  to  be  run  on 
a  CDC  6600  digital  computer.  It  is  capable  of  analyzing  configurations 
having  up  to  six  separate  "panels"  of  differing  flap  angle  and  jet 
momentum  per  side.  The  jet  distribution  and  orientation  can  be  varied 
so  as  to  correspond  either  to  internally  blown  flaps  (momentum 
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proportional  to  local  wing  chord,  blowing  normal  to  hinge  line)  or  to 
externally  blown  flaps  (momentum  unrelated  to  wing  chord,  blowing 
direction  parallel  to  the  aircraft  plane  of  symmetry). 

The  analysis  uses  a  nonplanar  trailing  vortex  system  in  order 
to  determine  the  component  of  induced  velocity  antiparallel  to  the 
relative  wind.  That  component,  neglected  in  the  analysis  methods  hitherto 
available,  becomes  important  when  the  very  large  circulation  lift  levels 
typical  of  jet  flaps  are  reached.  Its  effect  is  generally  to  reduce  the 
lift  corresponding  to  a  given  value  of  induced  drag,  so  its  neglect  is 
unconservative. 

Comparison  of  performance  data  predicted  by  this  analysis  with 
wind  tunnel  test  results  indicates  that  use  of  the  nonplanar  vortex 
system  concept  accounts,  in  several  cases,  for  substantial  drag  incre¬ 
ments  which  had  previously  defied  explanation.  However,  as  presently 
formulated,  the  procedure  over-predicts  drag  at  very  large  flap  deflec¬ 
tions.  Also,  some  difficulties  in  the  iteration  process  have  been 
encountered  where  local  jet  momentum  is  very  high  (i.e.,  in  externally 
blown  cases)  and  the  jet  deflection  angles  are  large. 

It  is  the  view  of  the  author  that  nonplanar  trailing  vortex 
geometry  is  an  essential  element  in  the  explanation  of  jet  flap  drag 
behavior,  and  that  the  method  of  this  report  deserves  further  develop¬ 
ment.  It  is  felt  that  adjustment  of  the  relation  between  jet  momentum, 
angle  of  attack  and  jet  angle  to  the  lifting  circulation  and  the 
trailing  vortex  geometry  details  can  lead  to  a  calculation  procedure 
giving  good  Results  for  the  whole  range  of  jet  flap  configurations 
of  interest  in  the  STOL  field. 
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SECTION  II 


THEORY 


2.1  Background 

The  theory  of  the  jet  flapped  wing  in  two  dimensions  is  complicated 
by  the  fact  that  it  is  a  "mixed"  boundary  value  problem.  The  customary 
requirement  of  classical  airfoil  theory,  zero  net  normal  flow  at  the  airfoil 
surface,  is  only  part  of  it.  Downstream  of  the  trailing  edge,  the  same 
condition  must  be  met  for  an  unknown  jet  shape.  In  addition,  the  loading 
there  must  just  suffice  to  turn  the  jet-  This  problem  was  solved  by  Spence^^ 
using  "classical"  mathematical  methods  and  by  Malavard®  using  the  rheoelectric 
analogy,  with  essentially  equivalent  results. 

The  jet  flapped  wing  of  finite  span  is  much  less  tractable.  Maskell 
and  Spence’s^  formulation  of  the  problem  could  only  be  solved  by  drastic 
simplification,  which  amounted  to  reducing  the  wing- jet  system  to  an 
elliptically  loaded  lifting  line.  Malavard^  obtained  some  interesting  solu¬ 
tions  to  the  three-dimensional  lifting  surface  problem  using  the  rheoelectric 
analogy,  but  the  inconvenience  and  special  equipment  required  make  that 
approach  unsuitable  for  design  purposes.  More  recently,  Lopez  and  Shen® 
and  Lissaman^  have  reported  relatively  convenient  methods  of  analyzing  jet 
flapped  wings  of  arbitrary  geometry. 

All  of  these  methods  use  linearized  boundary  conditions,  and  do  not 
account  for  the  streamwise  component  of  induced  velocity  due  to  the  vertical 
displacement  of  the  trailing  vortex  system,  as  diagrammed  in  Figure  1.  This 
streamwise  component  causes  a  reduction  in  lift,  but  no  saving  in  drag, 
for  any  given  value  of  the  "bound"  circulation  on  the  wing.  The  drag  polar 
can  be  significantly  affected  at  the  lift  coefficient  levels  generated  by 
the  jet  flap,  especially  in  the  case  of  partial  span  blowing.  Furthermore, 
this  effect  implies  an  upper  limit  to  aerodynamic  (or  "circulation")  lift 
which  is  determined  by  span,  as  opposed  to  section  characteristics. 

Helmbold!s  analysis®,  which  accounts  for  the  effect  of  roll-up  on  the  angle 
of  the  vortices  downstream,  gives  a  limiting  value  of  1.9  times  the  aspect 
ratio  for  circulation  lift.  Lift  and  drag  data  for  a  blown  wing  at  very 
high  momentum  coefficient,  reported  by  Lockwood,  Turner  and  Riebe^,  agree 
with  this  result. 

2.2  Strip  Analysis 

The  approach  adopted  here  is  to  divide  the  wing  into  strips  for 
which  the  section  lift  and  moment  coefficients  depend  only  on  the  local 
angle  of  attack  (o l)  ,  flap  deflection  (6p)  ,  and  momentum  coefficient  (Cj)  . 

This  sacrifices  some  realism,  especially  for  low  aspect  ratio  configura¬ 
tions,  but  greatly  simplifies  the  analysis  and  reduces  computation  time. 

2.2.1  Vortex  Structure 

Each  strip  has  an  associated  "horseshoe"  vortex,  as  diagrammed 
in  Figure  2.  Induced  velocities  are  calculated  at  control  points  located 
at  the  nominal  quarter-chords  of  the  centers  of  the  strips,  which  coincide 
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Planar  Vortex  System 


Three  Dimensional  Vortex  System 


Figure  1 :  Effect  of  Vortex  Wake  Vertical  Structure  on  Aerodynamic  Forces 


2:  Vortex 


with  the  centers  of  the  "bound”  segments  of  the  horseshoe.  Each  trailing  vortex 
is  divided  into  two  segments.  The  first  one  represents  the  "near  wake", 
the  location  of  which  is  dominated  by  the  section  geometry  and  local  blowing. 

It  extends  from  the  bound  vortex  line  to  a  point  extended  some  distance  beyond 
the  trailing  edge  of  the  flap.  The  extension  distance  (E)  was  determined 
empirically  to  give  good  drag  correlation  at  varying  c  ; 

J 


E/c  *  l.Z 


(l) 


The  second  segment  of  each  trailed  vortex  represents  the  "far  wake",  extend¬ 
ing  to  infinity  downstream  at  an  angle  (OC^  )  to  the  freestream  wind  direction  t 
determined  as  follows: 


1.  A  rough  estimate  of  the  local  circulation  lift  coefficient 
(c^  )  is  determined  from  the  local  Cj,  £p,  and  oc  (assuming 
zero  induced  velocity)  using  the  section  lift  function 
(discussed  later),  and  corrected  for  finite  span  using  the 
following  formula,  based  on  ordinary  lifting  line  theory: 


(no  down  wash ) 


A 

A+2 


(2) 


A  is  the  aspect  ratio  of  the  whole  wing,  and  is  assumed  to 
be  appropriate  for  this  correction. 

2.  The  wake  angle  is  inferred  from  an  approximation  to  the 
rolled-up  vortex  value  implied  by  Helmbold’s  highly 
loaded  wing  theory^: 


—  0.243  sClastx.  (c£c//3A^ 


(3) 


(A  derivation  is  given  in  Appendix  I.) 

2.2.2  Section  Characteristics 
"Jet"  Lift  vs  "Circulation"  Lift 

The  lift  on  a  jet  flapped  wing  includes  both  an  aerodynamic  (or  "circula¬ 
tion")  component  and  a  direct  jet  thrust  component.  To  determine  trailing 
vortex  strengths,  it  will  be  necessary  to  distinguish  between  them.  Further¬ 
more,  it  must  be  recognized  that  aerodynamic  forces  also  act  on  the  jet 
behind  the  wing,  and  thus  contribute  additionally  to  the  circulation.  Figure  3 
indicates  these  relations. 


The  total  section  lift  coefficient  is 


=  cjtCM  (ot  +  SF) 


(4) 
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Figure  3:  Section  Lift  Breakdown 


1 


The  aerodynamic  force  on  the  jet  must  be  just  enough  to  turn  the  jet  momentum 
vector  to  the  far  downstream  angle  £,  so 


-  <£ 


i7 


(5) 


Therefore,  the  "circulation"  lift  determining  the  vortex  strength  is 


Cj£c  ~  C4-  Cj  £  (6) 

The  drag  may  be  inferred  from  the  extra  streamwise  momentum  passing  through 
the  Trefftz  plane:* 


=  —Cj  oo-o.  e  (7) 


The  section  must,  therefore,  experience  an  aerodynamic  drag  coefficient 

Cdc  W  ~  C^ej><L  e  +  (8) 

This  force  appears  in  the  form  of  low  pressures  at  the  leading  edge,  which 
will  be  in  a  strong  local  upwash  region. 

In  the  two-dimensional  case,  there  can  be  no  trailing  vortices  and  £ 
must  be  zero.  Then,  all  the  lift  is  circulation  lift  (in  the  sense  that  it 
is  associated  with  bound  vorticity)  and  all  the  jet  thrust  is  recovered. 

The  two-dimensional  analyses  used  as  the  basis  for  the  lift  function  used 
here  compute  C#  c  on  t*le  t>asis  that  £.  is  zero.  If  it  is  not,  the 
bound  vorticity  on' the  jet  sheet  goes  down,  implying  a  reduction  in 
That  reduction  is  ignored  here. 

To  determine  Cjc>  £  must  be  specified  as  well  as  0C ,  S-p  and  Cj. 

In  the  present  analysis,  which  uses  an  iterative  procedure  to  establish  the 
circulation  of  each  horseshoe  vortex  element,  the  downwash  velocity  (w)  at 
the  lifting  line  (computed  for  the  previous  iteration)  is  assumed  to  be 
doubled  far  downstream.  Ignoring  the  streamwise  induced  velocity  component 
for  this  purpose, 


*In  the  present  analysis,  friction  or  "profile"  drag  is  neglected  because 
the  other  components  (jet  thrust  and  downwash-rotated  lift)  will  over¬ 
whelmingly  dwarf  it.  Note  also  that  the  section  coefficients  are  defined 
in  terms  of  the  direction  and  dynamic  pressure  (q)  of  the  local  relative 
wind.  They  must  be  resolved  and  adjusted  for  the  freestream  direction 
and  q  when  total  wing  forces  are  computed. 
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£ 


Xasn.  (2  -ur/v) 


(9) 


In  the  initial  determination  of  an  approximate  CjLc  to  define  a  trailing  vortex 
far  wake  angle,  £  was  assumed  to  be  zero.  The  subsequent  correction  of  c$ 
for  aspect  ratio  should  approximately  correct  for  this. 

Local  Relative  Wind  Adjustments 


The  application  of  section  data  to  three-dimensional  wings  is  best 
done  by  resolving  the  flow  into  components  normal  and  parallel  to  the  lifting 
line,  as  shown  in  Figure  4*  The  effect  of  taper,  which  implies  a  different 
resolution  at  each  chord  station,  will  be  neglected. 

Consider  a  typical  strip  (indicated  by  shading).  Its  area  will  be 
cext  i  y»  where  cext  is  the  local  streamwise  chord*  and  £y  the  width  of  the 
strip.  Let  the  blowing  thrust  applied  to  the  strip  be  denoted  by  <£J. 
Normally,  the  direction  of  the  blowing  thrust  will  either  be  streamwise  (as 
In  the  case  of  externally  blown  or  upper-surface  blown  flaps)  or  normal  to  a 
nozzle  line,  usually  the  flap  hinge  (as  in  the  case  of  internally  blown  jet 
flaps) . 


Now  construct  a  plane  normal  to  the  lifting  line  and  passing 
through  the  control  point  of  the  strip.  The  "normal  strip"  (indicated  by 
the  broken  lines)  of  width  &yf  and  length  cext  will  have  the  same  area  as 
the  streamwise  strip,  except  for  a  small  error  due  to  taper,  which  will  be 
neglected.  The  normal  component  of  jet  thrust  will  be 


ST 

WOGM/IL 


—  SJ  CJXL  A 


(10) 


for  streamwise  blowing,  and 


SJ 


A/OXMAL 


=  SJ 


(11) 


for  blowing  normal  to  the  hinge  line. 


Let  the  component  of  the  wind  vector  in  the  normal  plane  (including 
both  frees tr earn  and  induced  components)  be  denoted  by  V^,  Let  oc  and  be 
measured  in  the  normal  plane. 

The  momentum  coefficient  for  the  normal  strip  will  be 


C_  = 


SJvo/tMAt. 


(12) 


or  "Fowler  flap"  action. 
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including  extension, 


Component  of  Wind 
in  Plane  of  Cut  A-A 
(Normal  to  Lifting  Line) 

View  A-A 
(Rotated) 


Figure  4:  Relative  Wind  for  Strip  Analysis 


The  downstream  angle  of  the  jet,  used  to  determine  the  circulation  lift, 
must  be  revised  to  reflect  sweep  effects,  replacing  equation  (9)  above  by 
the  following: 


£.  =  fia^rC' [ 2 '&/ VCsCrtL.  ( A.+/3 )]  (13) 

where  A  and  /3  are  the  sweepback  and  sideslip  angles,  respectively.  The 
plus  sign  applies  to  the  left  wing,  the  minus  sign  to  the  right. 

The  circulation  lift  coefficient,  c ^  (as  defined  in  the  preceding 
subsection),  may  then  be  calculated  as  a  function  of  Cj,  a,  e,  and  c£/cext> 
using  the  section  lift  expression  given  in  the  following  subsection.  From  Cj  , 
the  bound  circulation  (y)  corresponding  to  the  strip’s  loading  may  then  be 
found : 


x  =  %  K  <*t/2  <14> 

The  moment  coefficient  (cm)  of  the  loading  on  the  normal  strip  about 
its  quarter  chord  can  also  be  extressed  as  a  function  of  Cj,  <Sy,  ot  and  Cp/c^xt_. 
This  contributes  a  couple,  acting  about  the  quarter  chord  of  the  normal  strip, 
and  parallel  to  the  lifting  line,  equal  to 

^A/  =  4  /°  Vf  cext  fy  '  (15) 

2.3  Lift  and  Moment  Functions 

Initially,  the  circulation  lift  function  used  in  the  present  analysis 
was  simply  an  approximation  to  the  curves  published  by  Malavard^.  It  was 
found  that  this  under predicted  measured  changes  in  lift  curve  slopes  and 
lift  due  to  jet  deflection.  Some  of  the  discrepancy  could  be  ascribed  to 
thickness  effects,  since  Malavard’s  electric  analogy  data  were  for  a  zero 
thickness  airfoil.  Wygnanski-^  suggests  that  entrainment  effects,  absent 
in  any  potential  flow  analysis,  would  also  amplify  lift. 

The  approach  taken  here  was  simply  to  adjust  the  coefficients  in 
the  section  lift  formula  to  give  good  agreement  with  the  data  of  Lpckwood, 
Turner,  and  Riebe^.  The  following  expression  is  used  for  the  total  lift: 

c£  =  C  +Ac.  (16) 

c-  xj 
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where 


C0  =  2tt[<X  +  (0.3Z  +  1.155  SA.  )  6  ]  (17> 

0  Cext 

is  the  lift  for  zero  jet  momentum,  and 

A  Cj  -  (2.1<b  SF  +/. 09Zot )  cj iM  +  (ol+  Q  Cj.  (18) 

J 

is  the  extra  lift  (both  aerodynamic  and  jet  thrust)  due  to  the  jet. 

The  moment  coefficient  (about  the  leading  edge)  is  determined  by  an 
approximation  to  Spence’s  thin  airfoil  results-*-*  2; 


where 


and 


in  which 


C  =  C 


Ac. 


*»LB  ML£.0 


(19) 


.r  Z  JT  a  ~F7  ,  -  2.33  F, 


LE..O 


ext  cex£ 


(20) 


(21) 


„  —1-204-  Cj-  , 

F  B/.25c.-^/.5f/-e  ) 


(22) 


The  moment  coefficient  referred  to  the  quarter  c^xt  point  is  then 

=  C«nLe  +  Cjl/4.  j  (23) 

Figure  5  shows  how  the  above  formulae  compare  to  the  relations  given  by 
Malavard  and  Spence. 

Section  "Model" 


Judgment  must  be  used  in  specifying  the  flap  deflection  angle  so 
that  the  lifting  effectiveness  of  the  actual  flap  design  to  be  analyzed  will 
be  fairly  represented  by  the  simple  c ^  vs  Sp  relation.  Figure  6  indicates 
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the  "model"  represented  by  that  function.  The  correct  fepr esentation  of 
multiple  segment  flaps  and  the  influence  of  jet  spreading  on  the  effective 
angle  of  the  jet  at  the  trailing  edge  are  beyond  the  scope  of  this  report. 


Figure  5:  Section  and  Cm 
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Direction 


Figure  6:  Section  Geometry  Conventions 


SECTION  III 


ANALYSIS  PROCEDURE 


3.1  Wing  Description 

Basic  Reference  Geometry 

The  procedure  presently  programmed  applies  to  trapezoidal 
wings  having  a  plane  of  symmetry,  although  these  restrictions  are  not 
essential  to  the  analysis.  Unsymmetrical  blowing  and  flap  deflection/ 
extension  are  permitted.  The  chord  line  at  the  center  of  the  wing  lies 
on  the  x-axis,  and  the  origin  is  at  its  25%  point.  Figure  7  shows  the  rela¬ 
tion  of  the  wing  to  the  coordinate  axes. 

The  wing  is  completely  defined  by  the  parameters  given  in  Table  I. 
All  dimensions  are  expressed  in  semi-spans,  so  the  area  of  the  wing  is 

sw  =  4/A  (24) 

Its  root  chord  is 

cr  =■  4/AU-hX)  (25) 


and  the  mean  aerodynamic  chord  (for  moment  reference  purposes)  is 

1.5  Sw 


The  x  and  z  coordinates  of  the  quarter  mean  chord  are  then 


(26) 


Xr,  = 


Cr  ~C 


c/A  (I  -  A) 


XoyfX  A 


(27) 


A-  ( /-  A) 


r1 


except  when  ^=1,  in  which  case 
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Figure  7:  Coordinate  Axes  and  Wing  Reference  Planes 


TABLE  I 


WING  DESCRIPTION  PARAMETERS 


PARAMETER  SYMBOL 


Aspect  Ratio 

A 

Sweepback  Angle 

A 

Taper  Ratio 

> 

Dihedral  Angle 

r 

Tip  Twist 

0 

Number  of  Panels 

n 

max 

Span  Station  of  Inboard 

End  of  nth  Panel 

yendn 

Flap  Chord  Ratio  (Based 
on  Nominal  Chord)  of 
nth  panei 

(V^n 

Flap  Deflections  Right 
and  Left  Wing,  nth 

Panel 

^R,n>  *FL, 

Flap  Extension  Ratios 
(Cext/c)  Right  and 

Left  Wing,  nth  Panel 

eFR,n’  eFL, 

Blowing  Type: 
a)  "Internal11 

Constant  cj 
over  panel 

(No  Symbol) 

or 

b)  "External" 

Constant  J 
per  unit  span 

- 

Blowing  Direction: 
a)  Normal  to  hinge 
line 

(No  Symbol) 

or 

b)  Streamwise. 
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Xr  =  0  .5  xta/n.  A 
v* 


0  .5  r 


(28) 


Several  unit  vectors  which  will  be  useful  later  are  defined  next: 
Nr  and  Nl  are  unit  vectors  normal  to  the  right  and  left  wing  reference 
planes,  given  by 


N(  j  =  Ox.  +  y  xUsn.  r  -h  k  ccxx  H 


(29) 


where  the  upper  sign  applies  to  the  right  wing^  and  i j  ,  k  are  unit  vectors 
in  the  x,  y,  and  z  directions,  respectively.  Br  and  BL  are  unit  vectors 
parallel  to  the  lifting  lines,  given  by 

B ^  j  =  —  x  A+y  ccr^Aco^T' . ±  k cca.  A  am*  P  (30) 

(The  sign  conventions  chosen  here  give  the  B  vectors  the  same  sense  as  the 
bound  circulation.)  As  discussed  previously,  section  characteristics  are_ 
considered  to  apply  in  the  plane  normal  to  the  wing  -  that  is,  normal  to  B. 

I^t  wil^  also  be  convenient  to  define  the  unit  vectors  perpendicular  to  both 
B  and_  N,  pointing  "downstream"  in  the  plane  normal  to  the  wing.  Since  B 
and  N  are  already  mutually  perpendicular. 


D 


( > 


(31) 


Trailing-Edge  Device  Geometry 

The  wing  is  divided  into  "panels"  corresponding  to  the  divisions 
of  the  trailing-edge  flaps  or  blowing  arrangements.  Up  to  six  panels  may 
be  used.  The  number  six  was  chosen  to  provide  the  capability  to  analyze 
a  four-engine  externally  blown  flap  configuration.  Flap  deflection  and 
extension  are  independently  specified,  according  to  the  conventions  stated 
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in  Section  II.  The  flap  chord  ratio,  referred  to  the  extended  chord,  for 
the  nth  panel  is 


(cF/c  t  )  =  fer/c)  / (32) 


Because  the  wing  is  analyzed  in  slices  0.04  semi-spans  wide,  the 
ends  of  the  panels  must  fall  on  span  stations  which  are  integral  multiples 
of  0.04.  The  nominal  area  of  the  nth  panel  is  denoted  by  Sn,  and  is 
calculated  according  to  the  procedure  given  in  Subsection  3.3  (Eq.  40). 

3.2  Flight  Condition 

The  additional  data  required  for  complete  statement  of  the  problem 
are  the  descriptions  of  the  relative  wind  and  the  jet  blowing.  Since  all 
calculations  are  nondimens ionalj_  the  freestream  velocity  is  taken  as 
unity.  The  freestream  vector,  U,  is  given  by 


U  *=  CCXL  CO-4/3 

-y  co-4.  susn.  ft  +  k  AMO.  «.  w 

where  <X^  and  are  the  wing  angle  of  attack  and  angle  of  sideslip. 

The  total  thrust  of  the  blowing  system  is  defined  by  the  wing 
thrust  coefficient,  Cjw,  referred  to  freestream  dynamic  pressure  and  the 
area  of  the  whole  wing.  The  distribution  of  the  thrust  among  the  panels 
is  defined  by  fg^  n  and  fgL  n.  These  are  weighting  factors  apportioning 
the  jet  thrust  among  the  panels;  e.g., 


^  J.  -n 


n.'n 


(34) 


Thus,  a  "panel  Cj"  can  be  calculated: 


(  ).n 


c 


Ju 


W 


Si 


(35) 


f'Tl 


where  Sn  is  the  area  of  the  wing  panel  corresponding  to  the  nth  spanwise 
flap  portion,  and  Sw  is  the  area  of  the  whole  wing.  Deflection  angles 
of  the  jets  are  taken  to  be  equal  to  those  of  the  corresponding  flaps. 

3.3  Determination  of  Strip  Parameters 

Geometry 

Twenty-five  equal-span  strips  are  used  on  each  wing.  The  points 
on  the  lifting  line  at  the  edges  of  the  strips  have  the  vector  coordinates 
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(36) 


_  + _ ■*-  Bn. _ 

2  5  003.  A  c£n.r 


where  i  is  the  "index"  of  the  point,  beginning  with  1  at  the  centerline,  and 
going  up  to  26  at  the  tip. 

The  control  point  for  the  i^  strip  is  thus  at 


(37) 


The  nominal  chord  of  the  i^h  strip,  measured  at  its  center,  is 


(38) 


The  corresponding  extended  chords  are 


'ext 


=  e 


C  )  ^ 


C  )f*l 


(39) 


where  n  is  determined  by  the  panel  on  which  the  strip  lies.  The  nominal 
area  of  the  it^1  strip  is 


v5<-  =0.0-4-  c 

% 

The  nominal  area  of  the  nth  panel  is  thus 


where  the  range  of  i  covers  the  strips  in  the  panel. 
Blowing 


(40) 


(41) 


The  thrust  coefficients  applicable  to  each  strip,  referred  to  free- 
stream  dynamic  pressure  and  to  the  extend ed  strip  area,  are 


<>,*- 


J(. )/n 


for  "internal"  blowing.  For  "external"  blowing 


(42) 
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(43) 


_  )t-n 

Sc  •  e 


SJ. 


where  mn  is  the  number  of  strips  in  the  panel. 

The  jet  behind  a  particular  strip  will  actually  induce  lift  on  the 
adjacent  strips,  and  the  absence  of  a  jet  behind  an  unblown  strip  will 
reduce  the  induced  lift  on  an  adjacent  blown  strip.  Furthermore,  large 
abrupt  changes  in  strip  characteristics  cause  rapid  spanwise  variations 
of  loading  which  slow  the  convergence  of  the  iteration  procedure.  It  has, 
therefore,  been  found  expedient  to  "smear"  the  edges  of  the  jet  by  trans¬ 
ferring  1/3  of  the  jet  momentum  on  any  strip  at  the  edge  of  a  panel  to  the 
adjacent  strip  on  the  next  panel. 


The  total  wind  at  the  it^1  control  point  on  the  indicated  side  is 


Vu.i  ~  U  +  Wn. 


W  will  be  available  from  computations  executed  in  a  previous  iteration, 
details  of  which  will  be  given  later. 

The  magnitude  of  the  wind  velocity  normal  to  the  lifting  line  is 


v,=  l^x 


The  angle  of  attack,  measured  in  the  normal  plane,  is 


•  N 


)  + 


The  thrust  coefficient  used  for  section  C ^  and  calculation  is  then 


for  blowing  normal  to  the  hinge  line,  or 


for  streamwise  blowing.  (Equations  46  and  47  follow  directly  from  equations 

10-12.) 
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to  be 


The  downstream  jet  angle  (£)  referred  to  in  Equation  13  is  taken 


(49) 


Next , 


and  are  computed  using  the  results  of  Equations 


45,  46  (or  47),  48  and  32  in  the  lift  and  moment  functions  of  Subsection  2.3. 


The  circulation  of  the  i*"*1  horseshoe  vortex  is  then 


y<  >,+  ~  ,  .•  ^ 


C  ).  ^ 


C  ),A~ 


'ext 


( ),* 


A)/2 


(50) 


3.5  Iteration  Procedure 


The  set  of  vortex  strengths  constituting  a  solution  is  found  by  succes¬ 
sive  approximation.  The  first  set  is  calculated  assuming  zero  induced  velocity 
at  all  control  points.  Substitution  of  each  new  set  of^'s  directly  into  the 
equations  for  induced  velocities  will  not ,  however,  lead  to  a  solution.  That 
procedure  is  unstable,  giving  a  rapidly  diverging  sequence  of  vortex  strengths. 
Instead,  the  iteration  process  is  "damped"  by  use  of  a  weighted  average  of  the 
previous  two  solutions  for  the  next  trial: 


= d< , 


0- J)  7c 


OLO  (  ), 


(51) 


The  damping  factors  (different  for  each  strip)  are  determined  by  the  process 
described  in  Appendix  III. 

In  cases  where  very  powerful  blowing  is  applied  over  narrow  portions 
of  the  span,  such  as  externally  blown  flaps,  the  abrupt  change  in  circulation 
between  adjacent  wing  panels  was  sometimes  found  to  drive  the  computed  local 
induced  velocities  and  flow  angles  to  values  where  the  nonlinearities  of  the 
analysis  would  prevent  convergence  of  the  iterations.  A  smoothing  procedure 
was  therefore  added.  This  process  leaves  the  computed  results  almost 
unaffected  where  convergence  is  obtained  anyway,  and  extends  the  range  of 
conditions  for  which  answers  can  be  obtained  to  substantially  higher  local  Cj 
values.  The  smoothing  procedure  is  also  given  in  Appendix  III. 

After  each  set  of  y' s  is  calculated,  it  is  checked  against  the 
previous  set.  When  the  maximum  difference  of  any  one  from  the  correspond¬ 
ing  OLD  is  less  than  one  percent  of  the  average  of  all  the  y's>  the  solution 
is  considered  satisfactory  and  iterations  cease.  Computation  of  final  output 
data  (overall  force  and  moment  coefficients,  load  distributions,  etc.)  is  then 
started. 

3.6  Induced  Velocities 

Induced  velocities  are  computed  from  aerodynamic  influence  coef¬ 
ficients  and  the  weighted  average  vortex  strengths: 
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O.k 


25-  _ 

^  [ Ynew  •  c ±tf.  ~*~X/ev  7l,  c ), -s/7 

'3  'r 


The  aerodynamic  influence  coefficients  are  simply  the  velocities  induced 
at  the  control  points  by  the  horseshoe  vortices  of  the  various  strips, 
computed  for  unit  circulation.  For  instance,  TRtL,i,j  is  the  velocity 
vector  at  the  ith  control  point  of  the  left  wing  due  to  the  jth  vortex 
on  the  right  wing,  for^R^j  =  1. 

Influence  Coefficient  Calculation 

The  T's  are  the  sums  of  contributions  of  straight  line  vortex  seg¬ 
ments,  computed  according  to  the  Biot-Savart  law.  (Appendix  III  gives 
.details.)  Figure  8  shows_ the  segment  arrangements  and  notation,  using 
Tr  l  i,j  ils  an  example.  P  and  Q  have  already  been  defined  by  Equations  36 
am  ^7!  Hr>t1  is  the  unit  vector  parallel  to  the  semi-infinite  segments  of 
the  far  wake  portion  of  the  vortices  trailing  from  the  n^  panel.  It  is 
given  by 

&«, ~  = z  04X1  (a»  -  K.~_  K/n) 

+  J coa  Mn/3  <53 

+  k 

where  oc ^  ^^is  computed  for  the  panel  by  the  procedure  given  in  Paragraph 
2.2.1  above.' 

The  corner  points  near  the  trailing  edge  are  found  from 


-5/ 


where 


CjSjy  [a7S  ^ J  & 

+  k  c-js.  ,[(%)*  *|7  ^  Sr 


The  quantity  E/c  in  Equation  56  is  the  effective  chord_ex tension 
due  to  blowing  given  by  Equation  1.  Large  abrupt  variations  in  K  and  L  at 
transition  points  between  panels  have  been  found  to  cause  anomalies  in  the 
computed  loadings  and  to  slow  the  convergence  process.  Therefore,  the  L 
at  the  outboard  end  of  a  panel  and  the  K  at  the  inboard  end  of  the  adjacent 
one  are  averaged,  and  the  common  value  used  for  both. 
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The  velocity  Tr^l^j  is  then  the  sum  of  the  contributions  of: 

1)  The  semi- inf inite  segment  from  infinity  downstream  to  L  .. 

_  _  3 

2)  The  segment  from  KR  ^  to 

3)  The  segment  from  PR  ^  to  PR>j+i« 

4)  The  segment  from  PR^+1  to  LRjj. 

5)  The  segment  from  LR  ^  to  infinity  downstream. 

If  the  control  point  in  question  were  on  the  right  wing,  TR  R  ^  j 
would  be  found  using  QR>i  instead  of  QL,i>  and  the  bound  vortex  segment ’ (tne 
3rd)  would  be  omitted,  since  its  contribution  must  vanish.  The  Tl^r’s  and  ’ 
TLjL's  are  found  using  the  same  procedure,  with  obvious  changes  of  subscript. 

3.7  Forces  and  Moments 

Vortex  Forces 


The  aerodynamic  forces  on  the  wing  and  jet,  which  give  rise  to  the 
vortex  system  in  the  first  place,  may  be  inferred  using  the  Kutta-Joukowski 
theorem: 


SF  =p(vxr)u 


(57) 


where ^F  is_  the  force  acting  on  a  segment  of  bound  vortex,  p  is  the  fluid 
density,  V  is  the  wind  vector  across  the  vortex  segment  (assumed  constant, 
since  the  segment  is  short),  P  is  the  circulation  vector  of  the  segment, 
and  iX  its  length.  Nondimensionalizing  forces  by  freestream  dynamic 
pressure  and  semispan  squared,  the  force  on  the  ith  segment  will  be 

F<),i  *  Z[ %,J.  5,)]^  (58) 


and  the  segment  length  will  be 

=  0.04  ^  A  P  (59) 


The  total  aerodynamic  force  acting  on  the  wing/ jet  system  is 

therefore 


(60) 


25 


Jet  Forces 


Part  of  the  direct  jet  forces  acting  on  the  wing  are  actually 
included  in  the  vortex  forces  acting  on  the  wing-jet  combination.  The  jet 
forces,  as  defined  here,  correspond  to  the  jet  momentum  far  downstream  of 
the  wing,  as  indicated  in  Figure  9.  The  nondimens ional  jet  force  increment 
will  be 


*  J(U  Cjst;>  epU/'n 


(61) 


The  vector  contribution  of  the  strip  jet  force  will  then  be 


for  blowing  normal  to  the  hinge  line,  or 

*^C),jL  =  [  C>0<L  A  [N(, 


for  streamwise  blowing. 


The  total  jet  force  will  be 

zs . 


J-TOT=  S?  (JR,.  +  JL,J 


(63) 


(64) 


A.  =  I 


Moments 


Each  strip  contributes  a  moment  increment  in  two  ways: 

1)  The  product  of  the  force  on  the  strip  and  its  distance  from  the 
moment  center. 

2)  The  couple  associated  with  the  section  Cm. 

The  i^  strip's  section  contributes 


M 


=  c 


C  ).A- 


/vn 


V 


N(),i 


>S  CAeF 
A  hl  J,™ 
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where  CmQ  j_  was  computed  at  the  same  time  as  Cjjcq  .  according  to  the  formulae 
given  in  Subection  2.3,  (The  Ms  T  s  are  referred  to  £ne  freestream  dynamic 
pressure  times  the  semispan  cubed.) 
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Figure  9:  Jet  Force  Determination 
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The  total  moment  about  the  origin  of  coordinates  is,  therefore, 


Mwr  =  <65) 


Force  and  Moment  Coefficients 


The  forces  and  moments  are  then  expressed  in  the  customary  coef¬ 
ficient  form.  First,  they  are  resolved  in  the  stability  axis  system.  Unit 
vectors  in  the  x,  y  and  z  stability  axis  directions  are  given  by 


/ 

Sz  =  Sx  X  Sy 

The  three  force  coefficients  are  then 

ct  -  -  rTor  •  s2/sw 


(66) 


(67) 


CD=  -  '  5X/SW  (68) 

Cy  =  Fnr  •  Sy/, Sw  (69) 

The  moments  must  be  transferred  to  the  quarter  m.a.c.  as  well  as 
resolved: 

C-sm  ~L^ror  ~h(A~  ^  ^  ^vj  ^70'1 

C#  =/ mtot  +(z  x-c/4  ■+  kzc/4)  x  rTOT]  *  sx/2  5  (71) 

~L^ror  Xc/4.  k-zz/J)  X  F-ror-J *  F>z/2-  5^  (72) 
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SECTION  IV 


EXAMPLES  OF  APPLICATION 


4.  General 


This  section  presents  applications  of  the  method  to  a  variety  of 
problems  for  which  wind  tunnel  test  data  is  available  for  comparison.  The 
configurations  range  from  the  somewhat  academic  case  of  a  rectangular  wing 
with  uniform  internal  trailing  edge  blowing  but  zero  flap  chord  to  the  very 
practical  case  of  a  swept  back  wing  with  double  slotted  flaps  and  external 
underwing  blowing. 

4.1  "Pure”  Internally  Blown  Jet  Flaps 

9 

Lockwood,  Turner,  and  Riebe  reported  results  of  tests  on  a  rec¬ 
tangular  wing  of  aspect  ratio  8.4  with  full  span  trailing  edge  blowing.  The 
semi-span  model  was  constructed  by  truncating  the  trailing  edge  of  a  symmet¬ 
rical  wing  and  installing  an  air  supply  tube,  as  shown  in  the  sketch  at  the  top 
of  Figure  10.  This  tube  fed  a  plenum  exhausting  through  a  nozzle  tangent  to 
the  tube.  The  jet  was  forced  to  follow  the  outside  surface  of  the  tube  by 
the  Coanda  effect  until  it  reached  a  sharp  corner,  the  location  of  which 
was  used  to  set  the  jet  deflection  angle.  The  "flap  chord"  was,  therefore, 
effectively  zero. 


This  model  was  tested  to  very  large  momentum  coefficients  (up  to 
56),  and  gave  results  corroborating  the  limitation  of  circulation  lift 
predicted  by  Helmbold®.  Figure  10  shows  measured  drag  polars  for  this  model 
at  Cj  =  7,  at  several  jet  angles.  The  test  results  are  given  both  as 
reported  in  Reference  9  and  as  adjusted  for  the  pressure  deficiency  on  the 
base  of  the  wedge  used  to  set  the  jet  deflection.*  "Classical"  jet  flap 
theory  gives,  for  elliptical  lift  distribution, 
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=  -  c 


Cu 


i r  Aw  +2  Cj 


(73) 


(This  neglects  the  contribution  of  parasite  drag,  which  is  negligible  in 
comparison  to  the  range  of  drag  values  generated  by  induced  and  jet  effects 
present  here.)  This  equation  agrees  very  well  with  the  adjusted  test  data 
for  low  jet  deflection  angles.  However,  at  j  =  56°  and  86°,  the  test 
results  break  sharply  away  from  the  classical  parabola.  The  fact  that  the 
planform  and  blowing  distribution  are  rectangluar,  rather  than  elliptical, 
does  not  explain  this  behavior,  since: 


1)  The  differences  are  much  too  large;  and 

2)  The  effect  of  the  planform  difference  (assuming  that  the  load¬ 
ing  retains  the  same  shape  as  the  lift  builds  up)  would  be  to 
produce  another  parabolic  polar,  corresponding  to  a  reduced 
aspect  ratio. 

*The  correction  was  determined  using  a  base  pressure  coefficient  equal  to  that 
required  to  give  Cv>  =  -Cj  at  zero  lift  and  =  0. 
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Figure  10:  Drag  of  a  Rectangular  Jet  Flapped  Wing  With  Full  Span  Blowing  (Ref.  9) 
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The  deviation  from  parabolic  form  at  high  is  similar  to  that 
given  by  Helmbold!s  theory  (already  cited),  in  which  the  anti-streamwise 
induced  velocity  due  to  the  tilt  of  the  trailing  vortices  reduces  the  lift 
at  a  given  circulation  value,  without  reducing  induced  drag.  This  effect 
should  be  even  more  striking  in  the  jet  flap  case,  where  blowing  at  high  c 
carries  the  vortex  system  steeply  down  behind  the  wing. 

The  method  of  this  report  (which  will  henceforth  be  referred  to 
as  the  3DV  method,  for  "3-Dimensional  Vortex")  includes  this  effect,  and 
as  a  result  is  successful  in  predicting  the  break-away  from  parabolic  shape, 
as  will  also  be  seen  on  Figure  10.  The  loss  in  lift  at  the  largest  angle 
tested  ( 6j  =*  110°)  is  prematurely  predicted,  indicating  that  further  refine¬ 
ment  of  the  analytical  relation  between  Cj,  £j ,  and  vortex  system  geometry 
is  required.  No  such  loss  could  have  been  predicted  by  a  theory  based  on  a 
planar  vortex  system,  however. 

Figure  11  shows  Cl  vs  O t  for  the  same  conditions.  Here,  the  need 
for  improvement  of  the  vortex  geometry  model  is  reflected  in  the  reduction 
of  computed  lift  curve  slope  at  & j  =  86°  to  zero  and  below,  which  the  test 
data  indicate  is  approached  only  at  =  110°  or  beyond.  Agreement  at  lower 
jet  angles  is  fair,  indicating  that  refinement  of  the  Cl  -  &j  -  CX  -  Cj 
relation  is  also  needed. 

4.2  Partial  Span  Jet  Flaps 

The  same  model  was  also  tested  by  Gainer-^,  with  blowing  over  the 
inboard  half  of  the  span,  at  Sj  =  57°.*  His  results  indicate  that  the  three- 
dimensional  trailing  vortex  structure  is  even  more  essential  in  the  analysis 
of  systems  using  partial  span  flaps.  Figure  12  shows  the  polar  obtained  at 
Cj  =  1.79.  In  this  case,  the  difference  between  the  test  results  and  the 
"classical"  jet  flap  polar  is  dramatic.  Here,  a  substantial  fraction  of  the 
added  drag  may  be  ascribed  to  the  short  span  of  the  heavily  loaded  inboard 
portion  of  the  wing.  This  effect  accounts,  however,  for  only  half  the 
measured  additional  drag  at  Cl  =  5.0.  The  3DV  method  predicts  nearly  the 
full  increment  indicated  by  the  adjusted  test  data. 

Figure  13  shows  the  calculated  spanwise  circulation  distribution 
for  the  Cl  =  6.15  condition.  The  lift  and  drag  of  a  planar-vortex  system 
having  this  distribution  of  bound  vorticity  were  computed  using  a  method 
based  on  Rauscher's  analysis^.  The  jet  forces  were  added  vector ially  as 
indicated  by  the  diagram  at  the  upper  right  of  the  figure  to  determine  the 
"classical"  drag  corrected  for  span  load  distribution  effects. 

The  failure  of  this  model  to  achieve  Cl's  above  5.0  is  attributed 
to  stall  of  the  outboard  portion  of  the  wing,  where,  of  course,  there  is 
not  jet  to  produce  low  trailing  edge  pressures  to  maintain  attachment  of  the 
upper  surface  flow.  The  3DV  computer  program  fails  to  predict  this  stall, 
of  course,  since  it  is  based  on  potential  flow  section  characteristics. 

*Gainer  includes  the  jet-deflection-setting  wedges  contribution  to  the 
chord,  giving  a  slightly  higher  reference  area  and  reducing  the  nominal 
aspect  ratio  from  8.4  to  8.3. 
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Test  Data:  I 

- As  Reported 

-  Adjusted  (Note:  Adjustment -0 
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Figure  1 1:  Lift  of  a  Rectangular  Jet  Flapped  Wing  With  Full  Span  Blowing  (Ref.  9) 
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*:  Drag  of  a  Rectangular  Jet  Flapped  Wing  With  Blowing  Over  the  Inboard  Half  Span 
(Ref.  11) 
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Figure  13:  Circulation  and  Lift  Distribution  for  Rectangular  Jet  Flapped  Wing  With  Blowing 
Over  the  Inboard  Half  Span 
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Figure  14  compares  measured  and  predicted  lift  curves.  Again, 
the  need  for  refinement  of  the  CL  -0£  -  -  Cj  relation  for  section  charac¬ 

teristics  is  indicated. 

4.3  Jet  Flap  AMST  Model 

A  model  representative  of  a  tactical  transport  equipped  with  internally 
blown  trailing  edge  flaps  was  tested  in  the  STOL  TAI  program,  and  the  results 
reported  in  Volume  IV  of  the  present  series  by  Monk,  Lee  and  Palmer^.  This 
model  was  more  realistic  than  the  somewhat  academic  device  discussed  in  the 
preceding  subsections,  in  that  the  wing  was  interrupted  by  a  cargo-transport 
fuselage,  the  jet  was  blown  over  a  75%  span  flap  of  substantial  chord,  and 
nacelles  were  installed. 

Figure  15  shows  the  measured  longitudinal  aerodynamic  characteristics 
of  this  model  for  two  jet  flap  angles  at  Cj  =  0.6. 

Lift,  drag,  and  moment  values  were  also  computed  using  the  3DV 
program.  Body  carry-over  effects  were  assumed  to  be  adequately  simulated  by 
an  unblown  flap  of  the  same  chord  and  deflection  as  the  jet  flap,  covering 
the  portion  of  span  between  the  inboard  ends  of  the  actual  flaps.  (The  topic 
of  body  carry-over  for  powerful  high  lift  systems  deserves  a  study  in  itself, 
and  was  not  investigated  in  this  study.)  The  flap  angle  used  for  computation 
was  5°  higher  than  the  model  geometric  value,  because  the  jet  was  nearly 
parallel  to  the  flap  upper  surface,  which  was  substantially  steeper  than  the 
centerline  reference. 

Agreement  is  good  at  =  40°.  However,  at  =  80°,  drag  and  nose- 
down  moment  are  underestimated,  and  the  lift  is  too  high.  The  drag  deficiency 
indicates  that  not  enough  3DV  effect  was  registered  at  $ F  =  80°.  In  the  pre¬ 
vious  case,  at  about  10  times  the  jet  strength,  the  3DV  effects  were  exaggerated 
at  high  angles.  The  need  for  refinement  of  the  relation  of  Cj  to  jet  geometry 
is  again  apparent. 

4.4  Comparison  with  the  EVD  (Lopez-Shen)  Method 

To  provide  a  comparison  of  the  3DV  method  with  recent  jet  flap 
analyses  using  more  sophisticated  techniques  than  the  classical  analysis, 
it  was  applied  to  a  case  already  analyzed  by  Lopez  and  Shen6.  This  was  a 
model  tested  by  the  Royal  Aircraft  Establishment  and  reported  by  Butler, 

Guyett,  and  Moyl^.  This  model  had  a  full  span  jet  flap  of  9%  chord,  but 
the  blowing  was  interrupted  at  the  center  by  a  fuselage.  The  wing  section 
was  an  NACA  4424  profile,  resulting  in  a  jet  angle  about  20°  steeper  than 
the  nominal  centerline  reference  angle  of  the  flap. 

Figure  16  shows  drag  polars  at  3  Cj's  for  30°  flap  deflection 
(about  50°  jet  deflection),  as  measured  by  the  RAE  and  as  predicted  by  the 
3DV  method.  Agreement  is  quite  good.  The  "EVD"  (elementary  vortex,  Douglas) 
method  applied  to  the  Cj  =  2.13  case  underpredicts  drag  substantially.  (Ref. 

6  did  not  give  EVD  drag  results  for  the  other  Cj  values.)  This  further 
substantiates  the  need  for  a  three-dimensional  vortex  representation. 
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Figure  14:  Lift  of  a  Rectangular  Jet  Flapped  Wing  With  inboard  Half  Span  Blowing  (Ref.  11) 
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Figure  1 5:  Longitudinal  Aerodynamic  Characteristics  of  a  Jet  Flapped  AMST  Model  at  Cj  =  0.6  ( Ref.  13) 


Figure  16:  Drag  Correlation  of  3D  V and  EVD  Methods,  RAE  Jet  Flap  Model , 
dF  =  30°  (8j  =  50°) 
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Figure  17  shows  lift  and  pitching  moment  for  the  same  test  condi¬ 
tions.  Both  methods  give  good  correlation  at  the  lower  Cj  values.  (In  fact, 
the  precision  of  the  3DV  data’s  agreement  at  Cj  =  0.4  is  fortuitous.  Since 
no  representation  of  airfoil  camber  is  included  in  the  CL  vs  a-  function,  oc's 
are  only  good  to  within  an  additive  constant.)  The  3DV  method  predicts  a 
reduction  in  lift  curve  slope  at  higher  angles  of  attack.  This  is  due  to  the 
reduction  of  dynamic  pressure  at  the  wing  caused  by  the  anti-streamwise 
component  of  induced  velocity.  This  prediction  is  borne  out  by  the  test 
results.  It  is  not  indicated  by  the  EVD  method. 

Moment  correlation  of  the  EVD  method  is  better.  However,  no  attention 
has  yet  been  devoted  to  "calibration"  of  the  Cm  -  Cj  -  $  Y  ~  04  relation  in  the 
3DV  program. 

4.5  Externally  Blown  Flap 

Farlett,  Greer,  Henderson  and  Car ter ^  report  the  results  of  testing 
a  swept-wing  cargo  transport  type  of  configuration  equipped  with  full  span 
triple-slotted  flaps  blown  by  underwing  nacelles.  Figure  18  shows  this  model 
In  the  "spread  engine"  arrangement.  (It  was  also  tested  with  the  engines  in 
dual  pods.) 


The  3DV  program  was  applied  to  this  case  to  determine  its  capability 
to  handle  problems  of  maximum  complexity.  For  analytical  purposes,  each  wing 
was  divided  into  five  panels.  The  center  panel  represented  the  inboard  wing 
plus  the  body  "carry-over"  region,  the  second  and  fourth  panels  represented 
the  portions  of  the  wing  blown  by  the  engines,  and  the  third  and  fifth  the 
remainder  of  the  wing. 

The  Cj  used  for  calculation  was  reduced  to  reflect  "turning  losses" 
as  measured  statically.  The  ratio  of  total  force  experienced  by  the  model 
to  thrust  developed  by  the  nacelles  alone  was  0.83  at  40°  flap  angle. 
Therefore,  Cj  =  1.55  was  used  for  the  3DV  analysis  where  the  test  was  run 
at  Cj  <=  1.87. 


Figure  19  shows  computed  and  measured  tail-off  longitudinal 
characteristics  for  this  case.  The  test  data  show  a  break  in  the  lift 
curve  at  OC  =  10°,  indicating  flow  separation  over  some  part  of  the  wing, 
most  likely  the  outboard  region,  judging  from  the  corresponding  drift  of 
the  pitching  moment  toward  more  positive  values. 

The  predicted  lift  curve  slope  (below  the  aforementioned  partial 
separation)  is  too  low,  although  the  absolute  lift  level  is  close  at  04=10°. 
The  pitching  moment  is  roughly  correct  at  low  CL,  but  does  not,  of  course, 
register  the  pitch  up.  The  drag  agreement  is  excellent  below  stall. 

Again,  refinement  of  the  lift  function  is  clearly  needed,  but  these  results 
were  considered  quite  encouraging. 

Test  data  were  also  given  for  Cj  =  3.74.  Computation  attempted 
for  this  case  failed  to  converge,  however.  In  this  case,  a  section  Cj 
nearly  equal  to  9  is  reached  behind  the  outboard  nacelle.  The  exact 
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Figure  17:  Lift  and  Moment  Correlation  of  3DV  and  EVD  Methods ,  RAE  Jet  Flap  Model  5c  =  30P  (bj  =  5CP) 
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Figure  18:  Externally  Blown  Flap  Wind  Tunnel  Model  (Ref.  15) 
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mechanism  by  which  the  rapid  local  variation  of  section  Cj  prevents  con¬ 
vergence  is  not  understood,  and  must  be  resolved  if  this  analysis  is  to 
reach  its  full  potential  usefulness. 

Engine-out  calculations  were  also  run.  Lateral  characteristics 
are  compared  to  test  data  in  Figure  20.  No  engine-out  tests  were  run  with 
the  vertical  tail  removed,  so  the  comparison  is  somewhat  ambiguous.  The 
predicted  rolling  moment  due  to  engine  failure  is  on  the  conservative  side 
until  stall  occurs.  This  implies  that  the  3DV  method  (as  presently  formulated) 
underestimates  the  spanwise  extent  of  the  influence  of  the  still- operating 
inboard  engine.  Additional  work  is  therefore  indicated. 

Because  the  model  had  a  vertical  tail,  comparison  of  yawing  moment 
and  sideforce  are  dubious.  However,  if  the  difference  in  sideforce  is 
attributed  entirely  to  the  fin,  then  the  yawing  moment  can  be  corrected 
as  indicated,  to  give  good  agreement.  (Admittedly,  the  scale  of  the  Cy  plot 
given  by  Reference  15  is  too  small  to  be  read  accurately,  but  the  correction 
should  be  good  to  about  0.01  in  C^.) 

4.6  Concluding  Remarks 

The  3DV  method  is  able  to  predict  the  drag  of  jet  flapped  wings 
at  high  blowing  coefficients  and  jet  angles  better  than  any  method  not 
containing  the  essential  feature  of  a  nonplanar  trailing  vortex  system. 

Refinement  of  the  section  C £  -Cl-  Cj  -  /p  relation  is  needed  to 
improve  lift  prediction.  The  usefulness  of  the  method  would  be  improved  by 
resolution  of  difficulties  in  the  solution  procedure  so  that  EBF  systems 
can  be  analyzed  at  higher  Cjfs  than  presently  possible. 

Further  study  of  the  spanwise  lift  distribution  due  to  locally 
concentrated  blowing  is  also  needed. 


43 


Figure  20:  Engine  Out  Lateral  Characteristics  of  Externally  Blown  Flap 
Configuration  (Ret  15) 
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APPENDIX  I 


DOWNWASH  ANGLE  OF  A  HIGHLY 
LOADED  WING 


The  vortex  system  trailed  behind  a  wing  will  be  inclined  to  the 
freestream  wind  vector  because  of  the  downwash  field  of  the  wing.  The  angle 
of  inclination  will  be  greatest  near  the  wing,  but  as  the  influence  of  the 
bound  vortex  diminishes  with  distance  downstream,  the  angle  approaches  the 
value  due  to  the  trailing  vortices  alone. 

Consider  an  elliptically  loaded  wing  of  span  b  and  midspan  circu¬ 
lation  PQ  in  a  stream  of  velocity  U  and  density  /O.  The  downwash  velocity 
w  induced  at  the  wing  by  the  trailing  vortex  system  will  be 

-20  =  n>/2t>  u-u 

Far  downstream,  the  induced  velocity  will  be  twice  as  great,  and  the  downwash 
(and  vortex  inclination)  angle  will  be 

£.  =  '  (r'o/Uh)  (1-2) 


The  velocity  w  will  be  inclined  to  U  by  at  least  this  angle,  so  the  wind  component 
at  the  wing  in  the  freestream  direction  will  be  reduced  to 


V  -  U  ~  £ 

or 

The  lift  on  the  wing  is 

L’T  Pvr* 

so 

2 

where  A  is  the  wing  aspect  ratio  (b  /area) . 


(1-3) 


(1-4) 


(1-5) 
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This  formula  implies  an  upper  limit  to  of 

C,  =  /J~  =  0.355 A  (i- 

*-um  o  y  3 

reached  when  r^/Ub  equals  ^ 2/3 .  However,  Lockwood Ts  data^  indicate  that  jet 
flapped  wings  can  attain  circulation  Cj/s  more  than  twice  this  "limiting" 
value. 


Helmbold®  pointed  out  that  this  problem  is  resolved  if  the  rolling 
up  of  the  trailing  vortex  system  is  considered.  Vortex  sheets  are  generally 
unstable,  and  tend  to  roll  up  into  pairs  of  more-or-less  concentrated  vortex 
cores.  Momentum  conservation  arguments  show  that  the  sheet  from  an 
elliptically  loaded  wing  should  roll. up  into  a  pair  of  concentrated  vortices 
of  strength  at  a  spacing  of  <Vb/4  from  each  other.  Each  one  induces  a 
downwash  velocity  on  the  other,  given  by 


r 

IxT  - 


(1-7) 


so  they  must  be  inclined  to  the  frees tr earn  at  the  angle 

s'  =  (2Q/rr*hU) 


(1-8) 


Helmbold  then  infers  a  relation  between  Cl  and  e1  through  an  argument 
based  on  energy  conservation  in  the  Trefftz  plane.  The  same  relation  can  be 
obtained  in  terms  of  conditions  at  the  wing  itself  by  assuming  that  the  down- 
wash  velocity  at  the  wing  is  still  the  value  given  by  Equation  1-1,  but  the 
inclination  is  reduced  from  £  to  £ ? . 

The  wind  at  the  bound  vortex  is  then 


or 


V'  -  U-  Hr  £  ' 

V' 


(1-9) 


(1-10) 


The  lift  coefficient  will  then  be 


TT  A 

2 


(i-n) 
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The  corresponding  limit  to  Cl  is 


CLum  =  rrz  A/3/3  =  1.8994 A  ' 


(1-12) 


This  value  agrees  well  with  Lockwood's  observed  maximum  circulation  Cl- 
The  "rolled  up"  vortex  inclination  angle,  £f,  is  therefore  used  as  the 
far  wake  inclination  angle  in  the  present  analysis. 


The  equation  relating  e'  and  Cl  is  awkward: 


rr  3/4 

4 


Ain.  €.'(/—  JUn.  2 c  ') 


(1-13) 


It  was  found  more  convenient  to  approximate  the  inverse  of  this  relation  by 

£./  =  0.243  J  (C  L  //.&994AJ  (1-14) 


This  expression  is  accurate  within  one  percent  over  the  whole  range  of  Cl's, 
from  0  to  1.8994A. 
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(48  Blank) 


APPENDIX  II 


VELOCITY  DUE  TO  A  VORTEX  SEGMENT  ' 


The  Biot-Savart  law  gives  the  incremental  velocity  due  to  an  element 
ds  of  a  vortex  filament  of  circulation  P,  at  a  point  whose  position  from  the 
element  is  r: 


d  V 


r 

4  it  |  r  I  5 


ds  xr 


(ii-D 


Integrated  over  a  straight  line  segment  of  constant  circulation,  this  gives 


V  I  = 


r 


4-t  r  h 


(  c-o-a  O,  ~ 


c.o- a 


OJ 


(II-2) 


where  h,  0]_  and  02  are  as  diagrammed  in  the  sketch  below,  drawn  in  the  plane 
containing  the  vortex  segment  and  the  point. 


The  direction  of  the  induced  velocity  would  be  "out  of  the  paper" 
for  the  circulation  sense  indicated  by  the  arrow. 


The  procedure  as  programmed,  therefore,  computes  the  induced 
velocity  vector,  for  unit  P,  using_the  vector  cordinates  of  the  segment 
ends  (P^  and  P2)  and  of  the  point  (P3)  ,  as  follows: 


1)  Define 
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(H-3) 


2)  The  unit  vector  parallel  to  v  is 

n  = 

UV  /  A  x  &! 


3)  The  cosines  are  given  by 


^  __  A  '  £> 
IAIIBI 


-and 


CO'vl  ^7 


A-C 


mzi 

4)  The  perpendicular  distance  is 

h  —  /  Bj  jusn  6,  — 


lA_x&i 

/A/ 


(II-4) 


(H-5) 


(II-6) 


5)  v  is  then  the  product  of  j  v  \  from  Equation  II-2  and  U 


In  the  case  of  a  semi-infinite  vortex,  P^ 
is  defined  as  a  unit  vector  parallel  to  the  vortex, 
Then  Steps  3  and  4  become: 

and 

and 

P3  are  as_before, 
replaces  P3  as  an 

but  A 
input 

3) 

The  cosines  are 

/Q  _  _ 

A  •  B 

CO<L  C7,  — 

and 

(Bt 

(II-7) 

CCrtL  &£  ~ 

/ 

4) 

The  perpendicular  distance 

is 

J)  =  1 A  x 

Bl 

(II-8) 
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APPENDIX  III 

DAMPING  FACTORS  AND  SMOOTHING 


1«  Damping  Factors 

The  vortex  strengths  are  determined  by  a  set  of  weakly  nonlinear 
equations.  The  nonlinearity  is  introduced  by  the  inclusion  of  streamwise 
induced  velocities  in  the  determination  of  angle  of  attack  and  in  the  rela¬ 
tion  between  circulation  and  lift  coefficient.  To  obtain  damping  factors 
which  would  result  in  economical  convergence  of  the  successive  approximations, 
consider  the  situation  when  nonlinearities  are  ignored.  The  equations  defin¬ 
ing  the  vortex  strengths  are  then  of  the  form 


and 


(HI-1) 


Xi  “  Ct  c  V/z 

A* 


(IH-2) 


where  C ^  and  C  are  (local)  partial  derivatives  of  section  C^* 

with  angle  of  attack  and  flap  deflection,  y  the  circulations,  Wjj  the  change 
in  OC  at  the  ith  station  due  to  unit  circulation  at  the  j  th  station,  and  Oiy  &y 
c  and  V  are  geometric  angle  of  attack,  flap  deflection,  chord,  and  wind  across 
the  lifting  line.  (All  of  the  latter  four  may  differ  at  various  stations,  but 
are  independent  of  then's.) 

Def ine 


*1/  =  i  CLotc  t111-3) 

Then  the  equations  for  the  y1 s  can  be  written 

x  +  ^  y,  =  (iii-4) 

where  includes  all  the  terms  not  containing  y' s .  To  solve  this  set  of 
equations  by  iteration,  an  initial  set  of^'s  is  assumed,  and  a  subsequent 
one  computed  using  Equation  III-4,  rearranged: 

=  ~  ^  y/oa>  <III_5) 


^''Circulation”  Cl  that  is. 
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The  next  set  of  trial  y  !s  will  then  be  given  by 


y,  =  d  ■  Yr  +  (/-  djYt.  cm- 

where  the  d-^'s  are  the  damping  factors  we  seek. 

The  Kj[ j  1  s  generally  follow  the  following  pattern: 

1)  For  i  =  j,  Kjj  is  a  large  positive  number. 

2)  For  i  =  j,  Kij  is  a  negative  number.  Unless  i  and  j  are  nearly  equal,  K 
is  a  small  negative  number. 

This  implies  that  the  solution  will  be  dominated  by  the  terms  for  which  i  =  j  , 
In  the  limit,  neglecting  all  other  terms,  that  would  imply 


l*Kii 


(III-7) 


For  what  value  of  dj_  will  this  result  be  obtained  for  *  starting  from  an 

arbitrary  set  of  y^^^'s?  Replacing  the  summation  in  III-5  by  the  single 
term^jK^i,  substituting  the  result  in  III-6  and  equating  to  III-7  gives 


_  /?. 


OLD 


/+K 


If  d ^  is  taken  to  be 


cL. 


J_ _ 

I  +  Kjji 


(III-8) 


(III-9) 


then  Equation  III-8  is  satisfied  identically,  and  the  coefficient  of^yi 

*  1  V  JjJJ 

vanishes. 


To  apply  this  to  the  actual  case,  Kii  is  interpreted  as 


K  =  _  -M 


OLO 


(III-10) 


This  is  calculated  by  computing  a  set  ofy's  based  on  zero  induced  velocities, 
then  finding  the  change  due  to  the  velocities  induced  by  a  very  small  incre¬ 
ment  in  y . 

In  practice,  the  d-^'s  so  computed  have  been  found  to  be  too  large, 
leading  to  divergence.  However,  values  in  the  range  of  0.8  to  0.9  of  the 
theoretical  dj/s  usually  work  very  well.  To  prevent  divergence  in  the  case 
of  anomalies  caused  by  possible  large  differences  between  the  computed 
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derivatives  and  those  correct  for  the  "near  solution"  environment,  the 
maximum  difference  between  XiOLD  anc*  X  iNEW  tracked.  If  this  "excess" 
increases  between  successive  iterations,  the  d-j/s  are  all  multiplied  by 
0.8  for  subsequent  trials. 

2.  Smoothing 


The  smoothing  procedure  begins  by  a  transformation  of  the  spanwise 
coordinate: 


_/ 


4^/n. 


III-ll 


Fourier  cosine  series  are  then  determine  for  the  circulation: 


/  (^)  =  ^  -Hw  ***■  ry/Z) 


III-12 


Each  side  is  done  independently,  which  implies  a  series  for  which  the  B's 
vanish  for  even  values  of  n.  The  result  is  that  the  smoothed  (y)  must  have 
zero  slope  on  the  airplane  centerline,  which  is  only  necessarily  true  for 
symmetrical  cases.  However,  the  slope  at  the  centerline  is  generally  small , 
even  when  symmetry  is  not  present,  in  the  flight  conditions  of  interest 
for  STOL.  No  serious  distortion  of  the  results  is  therefore  expected. 


The  coefficients  of  the  series  are  determined  by  numerical  integra¬ 
tion  of  the  expression  . 


OXl(^^/2)  dlij. 


III-13 


Smoothing  is  then  done  by  determining  a  new  circulation  function 
using  only  the  first  six  terms  of  the  series: 


_  "v 

y< 'smooth 


m  =  / 

(Odd.) 


III-14 


This  procedure  is  based  on  concepts  used  in  classical  lifting-line 
theory,  as  presented,  for  example,  by  Rauscher^.  The  transformation  and 
the  use  of  Fourier  series  are  "natural"  for  the  problem:  In  the  case  of 
elliptic  loading,  all  the  Fourier  coefficients  except  the  first  one  vanish. 
Furthermore,  when  the  wake  angle  is  low,  the  lift  is  determined  entirely  by 
the  first  coefficient: 


BA 


III-15 
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and  the  induced  drag  coefficient  is  given  by 


III-16 
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APPENDIX  IV 


COMPUTER  PROGRAM 


1.  General 

The  analysis  procedure  has  been  coded  for  automatic  computation  in 
FORTRAN  IV  language  for  the  CDC  6600  digital  computer, 

2.  Program  Structure 

The  program  consists  of  a  main  program,  FLAPZ2,  and  four  sub¬ 
routines,  VICTOR,  WASH,  CEEL  and  SMOOTH.  Communication  between  the  main 
program  and  the  subroutines  is  mostly  by  means  of  argument;  COMMON  is  used 
only  for  subroutine  SMOOTH.  Figure  21  is  a  block  diagram  of  the  main  program. 

2.1  Subroutines 

2.1.1  VICTOR 

Subroutine  VICTOR  is  a  package  of  three-dimensional  vector  operations, 
with  arguments  A,  B,  *C,  D  and  E.  Different  entry  points  are  provided  for  dif¬ 
ferent  operations,  as  tabulated  below: 


Entry 

Operation 

VPLUS 

C  =  A  +  B 

VMINE 

C  =  A  -  B 

VCROSS 

C  =  A  x  S 

VDOT 

D  =  A  •  B 

VMAG 

D  =  |A|  ,  : 

SCALM 

C  =  DA 

SCALD 

B  =  A/D 

2.1.2  WASH 

Subroutine  WASH  computes  the  velocity  vector  due  to  a  unit-strength 
straight  line  vortex  segment  of  finite  or  semi- inf inite  length.  Its  arguments 
are  P,  A,  Q  and  W. 

Entry  WASH  finds  the  velocity  W  at  Q  due  to  a  vortex  extending  from 
P  to  infinity,  in  a  direction  parallel  to  A. 

Entry  SEG  finds  W  at  Q  due  to  a  vortex  extending  from  P  to  A. 

2.1.3  CEEL 

Subroutine  CEEL  computes  ’’circulation  lift”  and  section  pitching 
moment  in  accordance  with  the  formulae  of  Section  2.3,  Equations  16-23. 
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Figure  21:  F LAPZ2  Program  Block  Diagram 


2.1.4 


SMOOTH 


Subroutine  SMOOTH  contains  the  smoothing  procedure  given  in 
Appendix  III.  The  argument,  G,  is  the  array  of  y  fs  (either  left  or  right 
wing)  to  be  smoothed.  The  array  of  smoothed  y's  is  transmitted  back  to 
the  main  program  by  the  same  argument. 

COMMON  contains  arrays  F  and  T  which  are,  respectively,  weighting 
factors  for  the  numerical  integration  used  in  computing  the  Fourier  coef¬ 
ficients,  and  the  values  of  the  transformed  coordinate,  y,  corresponding  to  the 
control  point  locations.  These  are  computed  in  the  main  program. 

3.0  Program  Listing 

The  FORTRAN  coding  of  the  program  is  listed  on  Pages  61  through  97. 
A  partial  table  of  correspondence  between  FORTRAN  variable  names  and  the 
quantities  referred  to  in  the  main  body  of  the  report  is  given  in  Table  II. 

4.  Input  Data 

4.1  Format 


A  complete  set  of  input  data  comprises  thirteen  punched  cards. 
Card  1  is  simply  70  columns  of  free-form  alphameric  data  which  is  printed 
as  a  heading  at  the  beginning  of  each  output  case. 


Cards  2  through  13  contain  numerical  data  in  n10-field,  seven-digit" 
floating  point  format.  That  is,  the  first  70  columns  of  each  card  are  divided 
into  10  equal  seven-column  "fields".  Each  field  must  either  contain  a  number 
with  a  decimal  point  (location  optional)  or  blanks  only.  Columns  71  and  72 
are  not  used.  Columns  73  through  80  are  available  for  card  identification. 


4.2  Data  Sequence 

Card  1:  Title 

Card  2:  Field  1 

Field  2 
Field  3 
Field  4 
Field  5 
Field  6 


Card  3:  Fields  1  to  5 


Aspect  Ratio 

Sweep  Angle  (degrees) 

Taper  Ratio 
Dihedral  (degrees) 

Twist  (degrees) 

Blowing  Direction  Number 
(Any  number  ^  0.  means  blowing  parallel 
to  x-axis.  A  number  <■  0.  means  blowing 
normal  to  hinge  line.) 

Span  stations  of  inboard  ends  of  flap 
panels  2  through  6.  If  less  than  6 
panels  are  used,  remaining  values  should 
be  1.0. 
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TABLE  II 


VARIABLE  NAME  CORRESPONDENCE  ' 


Wing  Descriptors: 


Derived  Wing  Data: 


Vortex  System  Data: 


Vortex  Strengths: 


Blowing  Data: 


FORTRAN  Name 

Text  Symbol 

AR 

A 

DIHDRL 

r 

SWEEP 

A 

TAPER 

A 

FB 

Yend 

CF 

cp/ c 

FFR,  FFL 

?F 

DFR,DFL 

^F 

TWIST 

e 

ESS 

Sw 

CR 

cr 

CBAR 

c 

XCBAR 

xc/4 

ZCBAR 

?c/4 

UNR,UNL 

N 

AB0NDR , AB0NDL 

B 

D0WNR , D0WNL 

5 

PR,  PL 

P 

QR>QL 

<5 

CN0M 

c 

CHDR , CHDL 

cext 

TEIR.TEIL 

K 

TE0R.TE0L 

L 

EHR,EHL 

H 

G0R.G0L 

yOLD  or 

-  /new 

GR,GL 

y 

CJW 

CJ 

Jw 

BL0R , BL0L 

fB 

CJPR, CJPL 

Cj 

CJR.CJL 

CJs 
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TABLE  II  (Continued) 


FORTRAN  Name 

Velocities  and  Influence  Coefficients: 

U 

TRR , TRL , TLR , TLL 
WR,NL 

Forces  and  Moments: 

F0R,F0L 

TF 

FJR,FJL 

TJF 

T0TM 

CLIFT 

CDRAG 

CSIDE 

CPITCH 

CR0LL 

CYAW 


Text  S' 


F 

fTOt 

J 

JXOT 

MTot 

cl 

CD 

Cv 
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Card 

4: 

Fields  1  to 

6 

Flap  chord  ratios.* 

Card 

5: 

Fields  1  to 

6 

Right  wing  flap  extension  ratios.* 

Card 

6: 

Fields  1  to 

6 

Left  wing  flap  extension  ratios.* 

Card 

7: 

Field  1 

Number  of  sideslip  angles  to  be  analyzed. 
(Max imum  9.0,  minimum  1.0.) 

Fields  2-10 

Sideslip  angles.  (At  least  one  value 
must  be  given.) 

Card 

8: 

Field  1 

Number  of  angles  of  attack  to  be 
analyzed.  (Maximum  9.0,  minimum  1.0.) 

Fields  2-10 

Angles  of  attack  (At  least  one  value 
must  be  given.) 

Card 

9: 

Fields  1  to 

6 

Right  wing  panel  blowing  factors.* 

Card 

10: 

Fields  1  to 

6 

Left  wing  panel  blowing  factors.* 

Note:  The  panel  blowing  factors  give  the  relative  blowing  thrust 

for  the  panels,  and  are  normalized  by  the  program  to  give 
a  sum  equal  to  1.0.  Therefore,  any  convenient  numbers 
having  the  correct  relative  magnitudes  will  work.  To 
indicate  constant  local  Cj  (as  for  "internal"  blowing) 
use  negative  numbers.  Otherwise,  constant  thrust  per 
unit  span  is  assumed. 


Card 

11: 

Fields  1 

to  6 

Right  wing  flap  deflection.* 

Card 

12: 

Fields  1 

to  6 

Left  wing  flap  deflection.* 

Card 

13: 

Field  1 

Number  of  Cj's  to  be  analyzed. 

(Minimum  1 . ,  maximum  8 . ) 

Fields  2- 

-9 

Cj  1  s 

Field  10 

Next  case  card  number.  (See  explanation 

below. ) 

The  following  case  (if  any)  may  be  defined  either  by  a  complete  set 
of  cards  or  by  a  partial  set,  beginning  with  the  card  indicated  by  the  "next 
case  card  number".  The  intent  of  this  arrangement  is  to  permit  analysis 
of  a  hierarchy  of  case  parameters,  with  Cj  the  most  easily  varied,  flap 
deflection  next,  blowing  distribution  third,  relative  wind  fourth,  and  wing 
description  the  most  trouble  to  vary. 


*For  Cards  4,  5,  6,  9,  10,  11  and  12,  if  less  than  six  flap  panels  are  used, 
the  extra  columns  are  ignored. 
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PROGRAM  FLaPZ2<  I  NPUT  .OUTPUT  .TAPE5=  INPUT  »TaPE6  =  'JUTPuT  ) 

DIMENSIONS.  ETC. 


COMMON  F  (  25 ) *  T  (  2  5  ) 

DIMFNSION 

UNR ( 3  ) 

♦ UNL ( 3) 

»  FER ( 6 ) 

. FFL ( 6 ) 

1 

R  LOR < 6 ) 

*  BLOL ( 6 ) 

»  DFR ( 6 ) 

»DFL( 6) 

2 

PR ( 26 . 3 ) 

»PL( 26.3) 

»  OR ( 25  »  3  ) 

*  OL ( 2  5  *  3 ) 

7 

ABONDR ( 3 ) 

, ABONDL ( 3 ) 

♦GRlTEt  3) 

.GLEFT ( 3 ) 

4 

FFSRt  25  ) 

»  FFSL ( 25) 

. DOWNR ( 3 ) 

. DOWNL ( 3) 

5 

CHDR ( 25 ) 

»  CHDL ( 25 ) 

*  DR ( 25 ) 

*DL ( 25 ) 

<S 

C JPR ( 6 ) 

,C JPL ( 6  ) 

*C  JR ( 2  5 ) 

» C UL ( 2 5  ) 

7 

TEIRI25.3) 

.TEILI25.3) 

.  AREX ( 3 ) 

, AL I  NX ( 3  ) 

8 

TFORl 25.3) 

.TEOl( 25.3) 

*  EHk( 25.3) 

. EHL( 25 .3  . 

Q 

TRR (25.25.3) 

»  TRL ( 25.25.3  ) 

.GR (25 ) 

»GL ( 25  ) 

X 

TLR ( 25  *  25 . 3 ) 

♦  TLL ( 25*25.3  ) 

« GOR ( 25 ) 

»GOL( 25 ) 

1 

OAMPR ( 25 ) 

.DAMPL (25 ) 

♦  W  R  (  2  5  »  3  ) 

♦  W  L  (  2  5 . 3  ) 

7 

D Y  NPR ( 2  5 ) 

♦DYNPL ( 25 ) 

.AJETR ( 25 ) 

.AJETLt  25  . 

3 

VELRI 25) 

♦  VELL (25) 

*  BETAR ( 25 ) 

.BcTAL ( 2  5  . 

4 

ALFARI 25) 

.ALFAL ( 25 ) 

. TW ( 25 ) . 

<3 

FORI  25.3) 

, FOL (25*3) 

*FJR( 25.3) 

.F JL( 25 .3  , 

6 

DLOADR<  25) 

*  DLOADL  <  25 ) 

.SLOADRJ 25 ) .SLOADLl 25) 

7 

UPLODR ( 25 ) 

. uPlODL ( 25 ) 

. C  MR ( 25 ) 

,CML(25I  > 

8 

SECCLRI 25) 

. SECCLL ( 25 ) 

♦  SECCDR ( 25 )  . SECCDL (25) 

DIMENSION 

W  1  ( 3) ,  W2 1 3) . 

W3 ( 3 ) ♦  W4<  3 ) ♦ 

W3(3)^ 

ROTORG ( 3 ) 

dimension 

TI TLF ( 10) 

.  X  I  N  (  1  O  ) 

♦  DV 1(3) 

. DV  2 ( 3) 

I B^T  A ( 9 ) 

. ALEF ( Q ) 

*CJ( 8) 

*U(  3) 

.TEMP ( 3  ) 

* 


2  T  FMP 1 (3) 


. TEMP 2 ( 3  ) 


*  CUE  ( 3 ) 
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.Pi  (3  ) 


.  P  2  (  3  ) 


3P i ( 3 ) 

,  P  4  1  S  1 

, AYE ( 3  1 

*  T  JM 1 3 ) 

,TJF( 3) 

4  T  F  (  3  ) 

,  TM ( 3  ) 

, TSM( 3  ) 

,FJET( 3) 

, 

5  SAX ( 3  ) 

,SAY( 3) 

»  SAZ ( 3  1 

, TOTE ( 3 1 

,TOTM<  3  ) 

6 FB ( 7  ) 

,CF<  6  ) 

, ND I V ( 7 ) 

*  LL ( 25 ) 

,CNOM( 25 ) 

7CFS( 25 ) 

♦  SPt  6  ) 

,BLOSf 6) 

LOGICAL 

CON,  SYM 

NJ.S  =  1 
N  A  S=  1 
N  B  S  =  1 

P  I  =3.141  59263 
DT0R=PI /]80. 

UNL  ( 1 )=0. 

U  NR ( 1  )  =0. 

CALL  DATE ( DATE  1 »  DATF2 ) 

NP  AGE  =0 
NpXT=l 
DG=. 00001 

c 

C  INPUTS 

C 

1  READ(5*101>  TITLE 

101  FORMAT ( 10A7 ) 

?  READ( 5*102)  AR*  SWEEP »  TAPfR,  DIHDRL,  TWIST,  BLOWDR 

iO?  FORMAT ( 1 0F7.0 I 
AR 19= 1 ,Q*AR 
DO  41  1=1,25 

41  TW<I)=TWIST*QTOR*( .04*I-«02) 

3  RFAD(5,102)  XIN 
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F  R ( 1 > =0. 


F8<7)=1. 

DO  51  J=  2 » 6 
M  FP(  J)=XIN<  J-l  ) 

4  R  FAD ( 5*102  )  Cp 

*  R  F  AD  ( 5  »  1  02  )  FFR 

FF  R  FAD (5*102)  FFL 

5  R  FAD (5*102)  XIN 

NR=IFIX( XlNt 1 ) ) 

DO  56  J= 1 *NB 

F6  BFTA( J ) =  X I N ( J  +  ]  ) 

READ<5*102>  XIN 
N A= I F I  X ( X I N ( 1  )  ) 

DO  57  J* 1 *NA 
F7  A  L  FF ( J ) =X I N ( J  +  l ) 

9  R  FAD ( 5*102 )  BLOR 

Q  CONTINUE 

I F ( NEXT , NF . 1 0 ) 00  TO  58 
DO  14  N=  1  ♦  6 

14  B  LOR ( N) =  BLOR ( N ) *SUMR 
F8  CONTINUE 

READ! 5*102) BLOL 

10  READ(5*102)  DFR 

11  READ(F,102)  DFL 

12  READ( 5* 102 )  XIN 

N  J=  I  F  I  X  (  X I  N  (  ]  )  ) 

DO  62  J=l*  NJ 

6?  CJt J)=  XINI J+l ) 


63 


c 


r  WING  GEOMETRY 

r 

SWP=5WEEP*DTOR 
SNSWs  SIN(SWP) 

C  55W=  COS(SWP) 

T  NSW=  TAN(SWP) 

TND I =  TAN(DIHDRL*DT0R ) 

DO  301  K=1 *3 

PR ( 1 *K) =0, 

301  PL(1»K)=0» 
SNDI=SIN(DIHDRL*DT0R) 

DO  302  1=1*25 

PR ( I +1  *  1 ) =  »04*TNS W# I 

PL (1+1 *1 1=  PR (  I  + 1  ♦  1  ) 

PR ( 1+1 ♦?  >  =  .04*1 
PL ( 1+1 *2 ) =  -.Oa*I 
PR ( 1+1*3)=  «04*TND 1  *  I 

302  PL ( 1+1 *3 ) =  PR ( 1+1 *3 ) 

E  SS=4./AR 

DO  303  1=1*  25 

DO  303  K= 1 ♦  3 

QR  (  I  ,K )= ( PR ( I » K ) +PR (  1  +  1 ,K)  ) /2. 
30?  QL(I.K)=(PL(I,K)+PL(l+l*<))/2. 

DO  21  N  =  1  *  25 

?1  T(N)  =  A  S I N ( OR ( N  *  2 ) ) 

F  (  1  >  =  .5*<T< 1 )+T( 2  )  > 


DO  22  N  =  2  *  24 
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22  F  (  N  )  *  ,5*(T(N+1)-T(N-1) ) 

F ( 25) =  .5*(T(25)-T(?4)  )  +. 5* ( P I /2 . -T ( 2 5 ) ) 

C 

C  NORMALS  AND  BOUND  VORTEX  VECTORS 

r 

UNR(2)=-5NDI 
UNL<2>=  SNDI 
DO  SOS  K.=  1 1  3 
GRITE(K)  =  PR  (  2  »  K  ) 

30k  GLFFT(K)  =-PL<2*<) 

UNR( 3 ) =5QRT ( 1 .-SNDI#*2) 

U  NL ( 3 ) =UNR ( 3) 

C 

C  SPANWISE  DIVISIONS 

C 

DO  306  L  =  2  f  6 

006  N  D I  V ( L ) =  I F I x ( ( FB ( L )  +.001)/. 04)  +1 

ND I  V  < 1 )  =  1 
NDIV( 7) =26 

CALL  VMAGIGLEFT.  DVT.  DV2»  G0M,  G0M2) 

FCTR  =  l./GBM 

CALL  SCALMIGLFFT*  DV1»  ABONdL.  FCTR*  DUM) 
CALL  SCALM(GRITE»  DV1.  ABONDRt  FCTR*  DUM) 

C  R  =  4./(AR*( 1 ,+TAPER )  ) 

C 

C  ASSIGN  PANEL  INDICES  TO  STATIONS 

r 

CRAR=CR**2*< 1 .+TAPER+TAPER**2 ) /< 1.5*ESS ) 
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I F { TAPER. NE.l . )G0  TO  307 
X  CBAR=. 5*TNSW 
ZCBAR=. 5*TNDI 
GO  TO  308 
30->  CONTINUE 

ARM= ( CR-CBAR ) / ( CR* ( 1 .-TAPER  )  ) 

XC8AR=ARM*TNSW 
ZCBAR=ARM*TNDI 
308  CONTINUE 

DO  310  Lx  1  »6 
I NDIO  =  NDIV(L) 

I NDIU  =  ND I  V ( L  +  l >-l 
00  310  I = I ND 10*  INDIU 
310  LL ( I ) =  L 

NP AN=LL ( 25) 

CALL  VCROSSI APONDR.  UNR,  DOwNR,  DUMA,  DUMC ) 
CALL  VCROSSI ABONDL ,  UNL ,  DOwNL,  DUMA,  DUMC) 
NP  A  NP 1  =  MP AN  + 1 

r 

C  FLAP  ANGLES,  CHORDS,  ETC. 

C 

DO  410  J= 1  *  25 
L  L J=  LL ( J ) 

DR(J)=  DFR( LLJ)*DTOR 
OL(J)=  DFL ( LLJ )*DTOR 
FFSL ( J ) =  FFL ( LL J ) 

CF< ( J)=CF<  LLJ) 


c  T  A= • O4*  J  -.02 
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FFSR(J)=FFR(LLJ) 


CN0M<J)=CR*(1 ,-ETA* (1  .-TAPER )  ) 

CHDR( J)*CNOM{ J ) *FFR ( LL J ) 

A10  CHDL( J)=CNOM( J)*FFL(LLJ) 
r  NOMINAL  PANEL  AREAS 

C 

DO  All  L  =  1 » NPAN 
SP ( L ) =0 • 

I  1 =ND I V ( L ) 

I  ?  =  NH>  I V  (  L  +  l  )-l 
IFII2.LT. II)  GO  TO  All 
DO  A12  1=11.12 

A  1 2  SP ( L  )  =SP ( L  )  +CNOM  C I ) / 2 5 • 

All  CONTINUE 
C 

C  BLOWING  DISTRIBUTION  -  NORMALISATION  OF  FACTORS 

C 

SUMB=0. 

DO  A21  L  =  1 » NPAN 

A21  SUMB=  SUMB+BLOL(L)+BLOR(L) 

ifisumb.eq.o, )auMB=1. 

DO  A 22  L  =  1  *  NPAN 
B  L0L ( L )  =  BLOL(L) /SUMP, 

B LOR ( L ) =  BLOR(L)/SUMB 
A  2  2  BLOS(L)=  2 . *BLOR ( L ) 

r 

c  overall  cj 

r 
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450  C  J3  =  CJ(NJS) 

r 

C  LOCAL  CJ  S  (PANELS  FIRST) 

r 

DO  431  L  =  1 ♦ NPAN 

IF (SP(L) .EQ.O.)  GO  TO  432 

C  JPL ( L ) =  BL0L(L)*CJ3*ESS/SP(L) 

CJPR(L)=  BLOR( L)*CJ3*ESS/SP ( L) 

GO  TO  431 
432  CJPL<L)=0, 

C JPR ( L ) =0. 

431  CONTINUE 

r 

C  CJ  S  AT  INDIVIDUAL  STATIONS  (REFERRED  TO  ACTUAL  CHORDS) 

C 

DO  435  1=1.  ?5 

L L I  =  LL ( I  ) 

ENSTA  =  NDIVILLI+l)-  NDIV(LLI) 

C JL ( I)=0, 

C  JR ( I  ) =0 . 

IFIENSTA.LT.l.)  go  TO  435 
IF (SUMB.LT.O. )  GO  TO  436 

C  JL (  I ) =  SP(LLI)*CJPL(LLI  ) / ( ENST A*.04*CHDL (  I  )  ) 

C  JR (  I ) =  SP ( LL I  ) *C JPR ( LL I  )/(ENSTA*.04*CHDR(  I )  ) 

GO  TO  435 

436  CJL( I )=CJPL(LLI ) / FFL ( LL I ) 

C  JR (  I  )  =  C JPR ( LL I  ) / E  FR ( LL I  ) 

43*  CONTINUE 
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r 

r  SMOOTH  CJ  variation  at  panel  ends 

c 

CJSPR  =  CJR ( 1  )/3. 

CJSPL=CJL( 1 1/3. 

CJRCl)*  CJR(1 1-CJSPR+  CJSPL#FFR( 1 ) /FFL( 1 ) 

C  JL (  1  )  =  CJL(1>-CJSPL+  CJSPR*FFL<1)/FFR< 1) 

IF  (NPAN.LT, 2)  GO  TO  470 
DO  460  N  =  2  » NPAN 
N  N  =  ND  I V  (  N  ) 

CJSPR  =  C JR ( NN ) / 3  » 

CJSPL  =  CJL(NN)/3. 

CJSPRl=  CJR ( NN- 1 ) / 3 • 

C  JSPL 1 =  C JL ( NN- 1 ) / 3 . 

C  JR ( NN ) =  CJR<NN)-CJSPR  +  C JSPRl*CHDR < NN- 1) /CHDR < NN ) 

C  JL  ( NN )  =  cJL  (  NN)  ~c  J^PL  +  COSPLi*CHDl-  ( ) /CHD1-  ( •'*■'* ) 

CJR INN- i )  =  CJR(NN-i)  -CJSPRl  +  C JSPR*CHDR ( NN ) /  CHDK(NN-l) 
460  CJL(NN-l)  =  CJL(NN-i)  -CJSPLl  +  C JSPL*CHDL ( NN ) /  CHDL(NN-l) 
470  CONTINUE 
r 

C  RELATIVE  WIND 

r 

SOO  BATER=BETA (NBS)*DTOR 
SNB=  SIN(BATFR) 

C  SB=  COS(BATER) 

T  NB=  TAN ( BATER ) 

ALFER=ALEF( NAS)*DTOR 


S  NA  =  SIN(ALFER) 
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CSA  =  COS(ALFER) 
CSSMB=COS(SWP-BATER ) 


CSSPB=COS( swp+bater ) 

U( 1 )=CSA*CSB 

U ( ? ) =-csa*snb 

IJ(5)  =  SNA 
WM  I  N=-.95*U  (  1  ) 
r 

r  sym.mptRy  t  EST 

( 

S  YM=  « TRUE  « 

I F { BATER  .NE .0  •  )  SYM=, FALSE. 

00  501  L  =  1  *  6 

I F  <  BLOR ( L ) .NE.BLOL ( L )  )  SYM= .FALSE. 

IF ( DFR ( L ) .NE.DFL ( L ) )  SYM=. FALSE. 

I F ( FFR ( L  )  .NE.FFL ( L )  )  SYM=. FALSE. 

A01  CONTINUE 
r 

r  TRAILING  EDGE  POINTS 

r 

DO  551  J=l,  25 
LLJ=  LL(J) 

CFLL=CF(LLJ)+  1 , 2*C JL ( J ) ** . ?5 

c  PL  =  CNOM  (  J  )  *(  .  75-CF(LLJ)+CELL*COS(DL(  J)  )  ) 

CBLR=CF ( LLJ )+  1 , 2* C JR ( J ) *  *  .  2  5 

FPR=CNOM ( J ) *( .75-CF( LLj )+CELR*COS( 0R( j ) ) ) 

T  E I L ( J . 1 )  =  PLIJ.l  )  +EPL 

T  EOL ( J  *  1 )  =  PL ( J+ 1  *  1 )  +EPL 
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T  E I R  (  J  ♦  1 )  =  PR( J.  1  )  +EPR 

TE0R(J*1)  =  PR ( J+l  *  1  )  +CPR 

T  E  I L  (  J  *  2  )  =  PL ( J  *  2  )  -EPL*TNB 

T  POL ( J ♦ 2 )  =  PLIJ+1*  2  )  -  EPL*TN0 

T  ElR ( J  *  2 )  =  PR ( J ♦ 2  )  -EPR*TNb 

TEOR  <  J »  2  )  =  PR  (  J+  1  *  2  )  -  E'PR*TNB 

F  PL  =  CHDL(J>*  CELL  *SI N ( DL ( J ) ) /FFL < LuJ ) 

F  PR  =  CHOR(J)*  CELR  >SIN I  DR ( J ) ) /FPR ( LL J ) 

TEILIJ.3)  =PL(J,3)-EPL 
T  FOL ( J  »  3 )  =PL(J+1»3)  -EPL 
T  E IR ( J  ♦  3  )  =PR(J*3)"EPR 
TEOR ( J  *  3 )  =PR(J+1»3!  -EPR 
CONTINUE 
r 

C  SMOOTH  TRAILING  FDGP  POINT  VARIATION  AT  PANEL  ENDS 

DO  510  K  =  1  ♦  3 

TEIR(l»K)=(TElR(l.K)+TFIL(l,K))/2. 

3  10  TEIL(1.K)=TEIR(1*K) 

DO  531  N  =  2 »  25 
NN=N 

DO  531  K=  1  ♦  3 

TFIR(NN,K)=  (TEIR(NN*K)+  TEOR  (NN-i  ,i<w)  ) /2. 

TEOR ( NN- 1 ♦ K ) =  T  F I R <  NN  *  K  > 

TEIL(NN,K]=  (  T  r  I  L  (  NN  «  K. )  +  T  FOL  (  NN- 1  ,M  )  /  2  . 

TEOL ( NN-l ,K ) =TPIL ( NN,K) 

SW  CONTINUF 

C 
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c 


trailing  vortfx  direction  cosines 


CALL  VCROS5 (  U»  ABONDL »  TEMP,  DUMA  *  DUMBC) 
CALL  VMAG  (TEMP,  DV1,  DV2»  WNORL,  WN0RL2) 

CALL  VDOT  (U,  'JNL,  DV1,  WNWCPE,  DtjMBC) 

ALLEFT  =  ASIN (WNWCPL/WNORL ) 

CALL  VCROSS  (  U»  ARONDR*  TEMP,  DUMA,  DUMC  ) 
CALL  VMAG  (TEMP,  DV1  »  DV2»  WNORR,  WNORR2) 

CALL  VDOT  ( U»  UNR,  DV1*  WNWCPR,  DUMBC) 
SPaNF*  AR/(AR+?.) 

A  LR I TE=  ASIN( WNWCPR/ WNORR ) 

DO  561  L=1«NPAN 
CFFL  =  CF ( L  )/FFL( L ) 

DLEFT=  DFL(L>*DTOR 
C  JLEFT  =  CJPL ( L ) /( WNORL2*FFL ( L  )  ) 

DR  I TF=  DFR(L)*DTOR 
C  JR  I TE=  CJPR(L) /(WNORR2*FFR(L)  ) 

CEFR=CF ( L ) /FFR( L ) 


CALL 

CEEL( ALLEFT, 

cjleft , 

DLFFT* 

cffl  , 

0«,  CLLEFT, 

DUMA ) 

CALL 

CEEHALRITE, 

CJRITE, 

DR  I  TE* 

CEER, 

0*»  CLRITE, 

DUMA) 

C  LLE  FT  =  CLLEFT*WNORL2/FFL( L ) 

CLR I TE=  CLRITE*WN0RR2/FFR( L ) 
r 

C  AT  THIS  POINT,  WE  HAVE  A  PANpL  CL  BASED  ON  NOMINAL  PANEL  AREA 

C  AND  FREE  STREAM  0.  NOW  FIND  A  WAKE  VORTEX  ANGLE  BASED  ON  THE 

C  WHOLE-WING  ASPECT  RATIO. 

C 

CLLEFT  =  CLLEFT*SPANF 

72 


CLRITE  =  CLRITF*SPANF 
1 F (CLRITE.GT.AR19 ) CLR I T E= . 9 9*AR 1 9 
I F(CLLFFT.GT.AR19 ) CLL EFT= . 99 *AR 1 9 
OFLTAR  =  ASlN(CLRITE/( 1.9*AR) )/4.l7 
DELTAL  =  ASIN(CLLEFT/(1.9*AR> J/4.17 
COSADL  «  COS(  ALFER-DELTAL  ) 

COSADR  =  COS( ALFER-DELTAR) 

A  L I  N  X  (  1 )  a  COSADL*CSB 
A  L I N  X  (  2  )  =  COS ADL*5NB 
A  L INX ( 3 ) =  SIN( ALFER-DELTAL) 

A  R  E  X ( 1 ) =  COSADR*CSB 
AR  EX ( 2 ) =  COSADR*SNB 
ARE  X ( 3 ) =  SIN( ALFER-DELTAR) 

I  1=  NDIV(L) 

12=  NDIV<L+1)-1 
I F (  I2.LT, II )  GO  TO  561 
DO  562  1=  Il»  I? 

DO  56?  K  =  l»3 

EHL ( I .K  )  =  ALINX  <  K  ) 

K6?  EHR ( I ,Kl  =  AREX(K) 

661  CONTINUE 
r 

C  COMPUTE  RIGHT  WING  INFLUENCE  COEFFICIENTS 

C 

DO  699  1=1.  25 
DO  60 1  K  =  1 ,3 
6^1  CUEtK)  =  QR(  I  ,K  ) 


DO  698  J= 1 .25 
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no  6O2  K  =  1  » 3 

P 1  (  K  )  =  TE 1 R  <  J »  < ) 

P?(0  cPR ( J  *K ) 

AYF(K)  =EHR <  J.K  ) 

P3<0  =PR ( J+l *K ) 

60?  P  4 ( K )  =T  EOR (J.K) 

CALL  WASH (Pi*  AYE.  CUE.  W1 ) 


CALL 

SEG(P1 . 

P?f 

CUE  . 

W2  ) 

CALL 

SEGIP3. 

ID 

CUF. 

W3  ) 

CALL  WASH ( P4.  AYE.  CUE.  WA  ) 

00  603  K=  1  » 3 

TRR(I.J.K)  =  W?(K>+  W3MO  +  WA  (  <  )  “  W 1  (  K ) 

PI  (<)=  T  EOL (J.K) 

P  7  (  K  }  =  PUJ+l.K) 

AYE ( K ) =EHL (  J.K  ) 

P  3  (  K  )  *PL  (  J  »  K) 

60^  PA  (  K )  =  TEI  L  (  J  »  K  ) 


CALL 

WASH ( P 1 » 

AYE  * 

CUE  * 

W1  ) 

CALL 

SEG 

(Pl» 

OsJ 

CL 

CUE » 

W2  ) 

CALL 

5EG 

<  P 2* 

P3  ♦ 

CUE* 

W3  ) 

CALL 

SEG 

( P  3 » 

P4  t 

CUE  ♦ 

W4  ) 

call 

WASH (PA, 

AYE  * 

CUE  » 

W5  ) 

DO 

60A 

K=1  * 

604  TLRU.J.K)  =  W2(<)+W3(K)  +  WA<K)  +W5(»s]  -  WHO 

AOs  CONTINUE 

IF(SYM)  GO  TO  699 
r 

C  COMPUTE  LEFT  WING  INFLUENCE  COEFFICIENTS 


c 


DO  605  X  =  1  *  3 

605  CUE ( X  )  =  QL< I.X) 

DO  699  J=l»  25 
DO  606  K=1 *3 

4YF(X)=FHR( J.  K) 
PI ( K ) =  TE I R ( J  *  <) 
P?(K) =  PR( J,  X) 

P  3  (  K  )  =  PR ( J  +  l *  K) 

606  P4(KI*  TFOR(J.  X ) 


CALL 

WASH( Pi  * 

AYE  ♦ 

CUE. 

wi  ) 

CALL 

SEG 

(pi. 

P2  * 

CUE  . 

W  ?  ) 

CALL 

S  E  G 

(  P2. 

P3  * 

CUF  » 

W  3  ) 

CALL 

5  EG 

(  P3. 

P4* 

CUE  , 

lV4  ) 

CALL 

W  ASH ( P4  * 

AYE  * 

CUE. 

w  5  ) 

no 

607 

K=  1 

*  3 

TRL ( I  *  J*  X)  =  W?(K)+  W3(X)  +W4(K)  +W5(X)  -W1(X) 

AYF(K)  =  FHL( J*  K ) 

PI (<)  =  TEOL( J»  X) 

P  2 <  K )  =  PL<  J+l .  X  > 

P3(K)  =  PL(J*  X) 

607  P4<X)  =  TE I L ( J  «  X) 


CALL 

WASHt  Pl  . 

AYE  . 

CUE* 

Wl  ) 

CALL 

SEG 

(Pl. 

P2. 

CUE* 

W  2  ) 

CALL 

SEG 

( P  3  * 

P  4  > 

CUE* 

W3  ) 

CALL 

WASH( P4. 

AYE  . 

CUE, 

W  4  ) 

DO 

608 

X  =  1 . 

3 

=  W2<K>  +  W  3  l  X  )  +  W4(X)  ~  W  1  (  K. ) 
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608  T LL ( I  ♦  J,  X) 


ftqo  CONTINUE 


C  FIRST  CUT  VORTFX  STRENGTHS  -  RIGHT  WING 

r 

OO  799  J=  1,  25 

CJCALL  =  CJR( J) /WN0RR2 
CFCALL  =  CFS ( J ) / FF  SR ( J ) 

IFIBLOWDR.GE.O.)  CJCALL  = C JCALL*CSSW 

CALL  CEELI ALR I TE »  CJCALL,  DR(J)*  CFCALL,  0.*  CLRR,  DCM  ) 
GRR=.5*CHDR  ( J)*CSSW*CLRR*WNORR 
GOR I J ) »GRR 
GRI J)=GRR 
r 

r  FIRST  CUT  VORTEX  STRENGTHS  -  LEFT  WING 

r 

CJCALL  =  CJL ( J ) /WN0RL2 
CFCALL  =  CFSI J) /FFSLI J) 

IFIBLOWDR.GE.O.)  CJC.ALL=  CJCALL*CSSW 

CALL  CEEL ( ALLEFT »  CJCALL,  DL(J).  CFCALL,  o.»  CLHL,  DCM ) 

G  LR=  .5*CHDL<  J)*CSSW*CLRL*Wf\|ORL 
G  L ( J ) =  GLR 
GOLIJ)  =  GLR 
7QO  CONTINUE 

CALL  SMOOTH ( GL ) 

CALL  SMOOTH(GR) 

CALL  SMOOTH (GOL) 

CALL  SMOOTH (GOR) 


I 0NCF=0 
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I  CON  =  0 


NTRY=0 

t  *5  0  CONTINUE 

r 

C  INITIALIZE  ITERATION  PARAMETERS 

C 

NTRIED  =  0. 

FXCESO=10000. 

CON=. FALSE. 

noi  continuf 

N  TPY=NTRY  +  1 

IF (NTRY.EQ.l )  GO  TO  802 
I F (NTRY.GT. 100 )  60  TO  85 
NTRIED  =NTRIED  +1 
I F (EXCESS. GT.O. 2 )  GO  TO  809 
IF(NTRIED.GT.30)  GO  TO  875 
80o  CONTINUE 
r 

C  AVERAGE  vortex  STRENGTHS 

r 

DO  802  J=l»  28 

GOR(J)  =  DAMPR(J)*  GR(J)  + ( 1 .“DAMPR ( J ) )*GOR( J ) 
IF(SYM)  GO  TO  802 

G OL ( J )  =  DAMPL(J)*  GL(J)  +(1.*DAMHL(J) )*GOL( J) 

802  CONTINUF 
r 

C  COMPUTE  INDUCED  VELOCITIES 


r 
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00 


80"* 


1=  1*  25 


00  804  K=  1 »  3 

WL< I ♦  K  )  =  0  • 

WR(I,K)  =  0* 

00  805  J  =  l»  25 

DO  805  K  =  1 »  3 

WR(I»<)=  WR (I »  K  )  +  GOR< Jj*TRR< I »J*K) 

IF (SYM)  GO  TO  806 

W  R  (  I  » K  )  =  WR(I»K)+  GOL< J)*TLR < I .J,K) 

W  L ( I ,  K )  =  WL(I,K)+  GOL ( J ) *T  LL ( I ♦ J  *  K ) 

W  L ( I »  K  )  =  WL ( I  *  K ) +  GOR( J)*TRL< I »J*K) 

GO  TO  805 
506  CONTINUE 

WR (!♦<)=  WR  (  I  »  K  )  +  GOR  (J)*T|_R<  I  *  J  ♦<  ) 

8  0  c  CONTINUE 

IF ( .NOT. SYM)  GO  TO  803 
WL ( I ,1 )=WR( 1,1) 

W L( I , 2) =~WR ( I ,  ?  ) 

W  L ( I ♦ 3 ) =  WR ( I . 3  ) 

8  0  ■a  CONTINUE 

DO  841  1=1.25 

IF(WR(I,1).lT.WMIN)WR(I»1)=WMIN 
I  F(  WL(  I  ,  l)  .LT.WMINJWUI  ,  l  )=WMIN 
841  CONTINUE 
r 

r  COMPUTE  DAMPING  FACTORS  (FIRST  PASS  ONLY) 

r 

IF (NTRY.NF.l )  GO  TO  880 
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CONTINUE7 


00  870  J=l.  25 
00  876  K  =  1  »  3 
CHECK=WR (  J»K  > 

JF(ARS(CHECK).LT.0.?5)  GO  TO  876 
CH ECK =. 25* ABS( CHECK) /CHECK 
876  TEMP(K)=CHECK+U(K) 

CALL  VCROSSI  TrMPi  ABONdR*  TFMP2*  DUMA*  DUMC) 
CALL  VMAG(TEMP2*  DVl*  DV  2  *  DNORR.  DNORR  2 ) 

CALL  VDOT ( T EMP *  UNR»  DVl *  DNWCP  *  DUMC) 

DAL  =  AS I N ( DNWC  P /DNORR ) 

CJCALL=  CJR ( J) /DN0RR2 
CFCALL=  CFS ( J ) / FFSR ( J ) 

CALL  CFELtDAL,  C-JCALL,  dR<U>.  CFCALc,  0.*  DCL » 
GRR=*5*CHDR  <J)*CSSW*DCL* DNORR 
DO  872  K=  1  *  3 
CHECK=WR ( J  *  K ) 

I F( ABS(CHECK) .LT.0.25 )  GO  TO  872 
CHECK=.25*ABS( CHECK ) /CHECK 
TEMP ( K ) =CHECK+U( K ) +nG*TRR( J,J*K) 

CALL  VCROSSI  TFMP  ♦  ABONDR*  T  r MP2  *  DUMA*  OUMC) 
CALL  VMAGITFMP2*  DV1*  DV 2 ♦  DNORR,  DN0RR2) 

CALL  VDOT (TEMP,  UNR*  DV1 *  DNWCP  ,  DUMC) 

DAL  =  ASINIDNWCP/DNORR) 

C  JCALL=  CJR( J)/DN0RR2 

CALL  CEEUDAL,  CJCALL,  D«<U>*  CFCALL,  0«»  DCL. 


DCM) 


DCM) 
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GRD  =  .  5*CHDR  (  J  >*CSSW*DCL*DNORR 
OGDG= (GRO-GRR ) /DG 
D  AMPR  (  J  )  =  l./U.-DGDG) 

DAMPRI J)=ABS(DAMPR( J) ) 

DO  877  K  =  1 »  8 
CHECK=WL  <  J*K) 

I F ( ABSf CHECK) ,LT. 0.25 )  GO  TO  877 
CHECK-. 2 5*ABS( CHECK ) /CHECK 
877  TEMP ( K) =CHECK+U( K ) 

CALL  VCROSS ( TEMP ,  ABONDL,  TEMP2,  DUMA*  DUMC ) 

CALL  VMAGtTEMP?*  DV  1 ,  DV2  *  DNORL*  DN0RL2 ) 

CALL  VDOT ( TEMP  *  UNL,  DV1*  DNWCP*  C JMC ) 

D  AL=  A 5. I  N(  DNWCP/DNORL  ) 

CJCALL-  C JL ( J ) /DNORL 2 
CECALL=  CFS(J) /FFSL(J) 

CALL  CEEL  (DAL,  CJCALL,  DL ( J )  ,  C  FC ALL  *  0.*  DCL*  DCM ) 

GLR=  ,5*  CHDL ( J ) *CSSW*DCL* DNORL 
DO  877  K  =  1  *  3 
CHECK-WL (  J  ,K ) 

I F( ABS(CHECK) .LT.0.25 )  GO  TO  873 
CHECK  =  ,2  5*ABS<  CHECK) /CHECK 
873  TFMP(K)=CHECK+U(K)+DG*TLL( J,J*K) 

CALL  VCROSS ( TEMP  *  ABONDL*  TEMP2»  DUMA,  DUMC) 

CALL  VMAG(TEMP2,  DV 1  *  DV2 ,  DNORL*  DN0RL2 ) 

CALL  VDOT (TEMP,  UNL*  DV1*  DNWCP*  DUMC) 

DAL=  AS  I N ( DNWrP/ DNORL ) 

CJCALL=  CJL ( J5/DNORL2 
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CALL  CEFL  (DAL,  CJCALL,  DUJ), 


CFCALL , 


0 . ,  DCL , 


DCM  ) 


GLD=  .5*  CHDL ( J)*CSSW*DCL*DNORL 
DGDG=  <  GLD-GLR ) /DG 
D  A  MPL ( J )  =  1 • / ( 1 » -DGDG ) 

DAMPLt J)=ABS(DAMPL( J) ) 

87o  CONTINUE 

OAMREF=DAMPR< 1 ) 

880  CONTINUE 
r 

r  COMPUTE  NEW  VORTEX  STRENGTHS 

C 

00  810  1=  1,  2  8 
DO  811  K=  1 , 3 

«11  TEMP(K)  =  WR ( I » < )  +  U(K) 

CALL  VCROSS ( TEMP  *  ARONDR,  TEMP2*  DUMA,  DUMC ) 

CALL  VMaG(TEMP2»  DV 1 ,  DV2,  vNGRR,  v,'iORR2) 

CALL  VOO T ( TEMP ,  UNR,  DVl ,  WNWCP,  DUMC) 

DYNPR ( I ) =VNORR  2 

A  LF=  ASIN ( WNWCP/VNORR)  +TW(I> 

C JCALL=  CJR ( I  I /VN0RR2 
I F (BLOWDR.GE.O. )  CJCALL=  cucall*cssw 
ALFIN=ATAN( -2.*WR ( 1,3) /CSSMfl ) 

CFCALL  =  CFS<  I ) / F  F  SR (  I ) 

CALL  CEEL ( ALF ,  CJCALL,  DR ( I ) ,  CFCALL,  ALFIN,  CLRR , 
A JFTR ( I ) =  ALFI N-ALF 

G R  t I )  =  . 5*CHQR (  I  )*CSSW*CLRR*VNORR 

C  MR ( I )  =  CMR (  I  ) +  CLRR*( l.-FFSRl I )  ) /4, 

IF (SYM)GO  TO  810 


CMR(  I  )  ) 


DO  812  K=  1  *  3 
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8  1  ? 


TFMP(K)=  WL(T  »<)  +U(K  ) 


CALL  VCROSStT^MP*  ABONQL,  TFMP2»  DUMA*  DUMC) 

CALL  VM  AG<  TFMP2  *  DV1*  DV2*  VNORL,  VNORL2 > 

CALL  VnCT(TFMP,  UNL  *  DV  1  ,  WNWCP*  DUMC ) 

DYNPL(I)  =VNORL2 

ALP=  A5IN( WNWCP/VNORL)  +TW(I) 

CJCALL  =CJL 1  I  ) /VNORL2 

I  FOLOWOR.GF.O.)  CJCALL  =CJCALL*CS$W 
ALFIN  =  ajAN(-2.*WL(  I  *  3 ) /C SSP 8 ) 

A JETL  < I )=ALFIN-ALF 
CFCALL  =  CFS (  I  > / F  F  5  L (  I  ) 

CALL  CEELt  ALF*  CJCALL*  DL(I)»  CFCALL,  ALFIN,  CLRL.  O'iL  (  J  ,  ) 
CML(I)=  CML ( I )  +  CLRL*( l.-FFSL ( I ) ) /4. 

GL ( I )  =  .5*  CHDL ( I )*CSSW*CLRL*VNORL 
310  CONTINUE 

TF(SYM)  GO  TO  815 
CALL  SMOOTH (GL) 
a18  CALL  SMOOTH (GR) 
r 

r  COMPUTF  EXCESS 

r 

E  X  C  E  S  5  =  0. 

SUMG  =0. 

DO  820  1=  1*25 

gl'MG  =  SUMG  +  GR  (  I  ) 

D I F  =  ARS ( GR (  I )-  GORI  I  )  ) 

I FIDIF.GT. EXCESS)  FXCESS  =DIF 


I F ( SYM )  GO  TO  820 
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SI  IMG 


SUMG  +  GUI) 


I  F  =AR.S(6U  I  )-  GOL  (  I  )  ) 
IF(DIF.GT. EXCESS)  EXCESS  =  DIF 
p?0  CONTINUE 

GAVE  =  SUMG/50. 

IF (SYM)  GAVE  =  SUMG/ 2  5 • 

EXCESS  =  EXCESS/GAVF 
EXCESS=ABS<  EXCFSS) 

IF (EXCESS.LT. EXCF SO)  GO  TO  83 
I  CON=  0 

DO  8 3 T  1=1*25 
OAMPR ( I ) =,8*DAMPR (  I  ) 

o'M  0 AMPL (I )  =  .8*DAMPL (  I  ) 

RATI 0=DAMPR ( I ) /  DAMREF 
GO  TO  750 
«30  CONTINUE 

CONTINUE 
EXCES0=EXCESS 

IFIEXCESS.GT..01 >  GO  TO  801 
860  CON=  .TRUE. 

880  NTRY=NTRY-1 

C 

C  COMPUTE  VORTEX  FORCES 

C 

00  901  1=1.25 

GR ( I ) =  «  5* ( GR ( I ) +  GOR (  I  ) ) 

DO  002  K  =  1  »  3 

oOp  TEMPI K)  =  U(K)+  WR(I.K) 


CALL  VMaG(TEWP»DV1 *  DV2  *  VELR(I),  DUMC ) 

CALL  VCROSS ( TEMP  *  GR I TE »  TEMPI  *  DUMA*  DUMC) 

DO  9010  K  =  1  *  3 

A  L  F  AR ( I )  =  AS  I  N  (  TEMP  (  3  )  /  R  T  (  TEMP  (  i  )  **2+TEMP  (31**2)  )  /DT^R 
ALFARt I )=ALFAR( I )+TW( I )  /DTOR 
oOlO  F  OR ( 1  *  K)  =  TEMPI ( K )  *  GR(I)  *  2. 

«n1  RFTAR(I)=-  ASIN( TFMP(2)/VELR( 1 ) (/DTOR 
IF (SYM)  GO  TO  910 
DO  910  1=1,25 
G  L ( I  )  =  ( GL (  I ) +GOL (  I )  > / 2 • 

DO  004  K=l,  3 
004  TEMP  (  K. )  =  U  (  K  )  +  WL(I,K) 

CALL  VCROSS ( TEMP  *  GLEFT,  TEMP2*  DUMA*  DUMC) 

CALL  VMAGITFMP,  DV1*  DV2  *  VeLL(I),  DUMC) 

DO  903  K=  1  *  3 

F OL ( I  ♦  K )  =  TEMP2  (  K  )  ■*  GL<I)*2. 

onu  CONTINUE 

PFTAL ( I ) =  -  ASIN ( TFMP { 2 ) /VELL ( I ) ) /nTOR 

4LFAL ( I ) =  AS  I N ( TFMP ( 3 ) /SORT ( TEMP ( 1 ) **2+TEMP  ( 3 ) **2 )  J/DTOR 
ALFAL ( I) =ALFAL ( T ) +TW( I )  /DTOR 

oio  CONTINUE 
r 

C  SUM  FORCES  AND  MOMFNTS 

r 

DO  920  K= 1  *  3 
T  F  {  K  )  =0  • 

TM(K) =0* 
opn  T  SM ( K ) =  0 • 
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DO  921  1=1.  25 

DO  921  K= 1  ♦  3 
o?i  TF(K)  =  TF  (  K. )  +  FOR  (  I  *  K  ) 

I F ( , NOT , SYM )  GO  TO  930 

CALL  SCALMITF.  DV1»  TF.  2.»  DUMC) 

DO  922  1=1*  25 

*2?  TV ( 2 ) =TM ( 2 ) +FOR ( I , 1 )*QR( 1*3) -FOR( I  * 3  > *QR (1*1) 
TV(2)=2.»TM(2) 

GO  TO  94 
<?30  CONTINUE 

DO  931  1=1*  25 

DO  932  K  =  l*  3 
T EMP  (  K. )  =  QR  <  I  *  K  ) 
op  TFVPl  (  K  )  =  FOR  (  I  , K  ) 

CALL  VCROSS(TFMP*TEMPl,  TEMP2.  DUMA*  DUMC) 
CALL  VPLUS ( TEMP2  *  TM,  TM,  DUMA.  DUMC) 

DO  933  K=1  .3 

TEMP(K)  =  OL(I.K) 

933  TEMPI ( K I =  FOL( I .<) 

CALL  VCROSS ( TEMP «  TEMPI,  TFMP2,  DUMA.  DUMC) 
CALL  VPLUS ( TEMP 2 .  TM,  TM.  DUMA.  DUMC) 

DO  934  K=  1  »  3 

934  TF(K)=TF(K)+FOL( I *K) 

931  CONTINUE 

040  CONTINUE 

r 

r  COMPUTE  jft  FORCES 

r 
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IPIBLOWDR.LT.O.)  GO  TO  945 
DO  941  1=1 *25 

5TR I PJ=CHDR ( I )*CJR( I) *.04 
DO  942  K=  1  *  5 

FACTOR=  -ABONDRf K )*SNSW  +  CSSW* ( S I N < A JE TK ( I ) ) *UNR { K ) 

1  -COS( AjETR ( I  )  )*  DOWNR ( K )  ) 

Q42  FJR(1»K)  =  $TRi PJ*FACTOR 
041  CONTINUE 

I  c  (  S  YM  )  GO  TO  950 
00  943  1=1*  ?5 

S  TR  I  P  J=  CHDLtl)*  C  JL  (  I  )  *  «  04 
DO  944  K=l*  3 

FACTOR=  ABONDLt K)*SNSW  +  C SSW* ( S I N ( AOETL ( I ) j *UNL( K ) 

1  -COSt AJETL ( I ) )*  DOWNL ( K ) ) 
o  44  F  J  L  (  I  *  K )  =  STR I P  J*  FACTOR 
P4^>  CONTINUE 

GO  TO  950 
QuS  CONTINUE 

DO  946  1=1*  25 

S  TR I P J=  CHDRt  I)*CJR(  I)*«04 
DO  946  K  =  1  *  3 

045  cjR(ItK)=  STRTPJ*(SIN(AJETR(  m*UNK(K.)  -  COS  (  A-l  ETR  U  )  )  *DOWimK  (  K  )  ) 
I^tSYM)  GO  TO  950 
DO  947  1=1*  25 

STR I PJ=  CHDL ( I ) *  CJL(I)*,04 
DO  94T  K=l*  3 

ozjT  cJl(j,K)-  STRtPJ*(S!N(aJETL(  I  )  )*UNl(M  -  COSt  AJETLt  I  n^DOWiMi-tN)  ) 
050  CONTINUE 
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c 

C  SUM  jet  forces  and  moments 

c 

00  949  K  =  1  *  3 

T JF ( K) =0. 

°4°  TJM(K)  =  f). 

DO  951  1=1*  25 

DO  952  K  =  1 »  3 

TEMP ( K ) =  QR (  I «  K ) 

959  F JET  <  K) =  F  JR (  I ,K ) 

CALL  VCROSS(TEMP*  FjET  »  T  cMP 1 »DUma ,  DuMC) 
CALL  VPlUS(T-»F*  FJET,  TJF,  DUMA*  DUMC) 
CALL  VPLUSITJM,  TEMPI,  TJM,  DUMA,  DUMC  ) 
o«?l  CONTINUE 

IF  (  S  YM  )  GO  TO  959 
DO  955  ! =1 ,25 

DO  954  K=1 *  A 

TCMP(K)  =QL  (  I  »  K. ) 
o  R4  FJFT(K)  =FJLU»K) 


CALL. 

VcROSSt  TEMP  , 

FJET, 

TEMPI 

,  DUMA* 

DUMC  ) 

CALL 

VPLU5 

(  T  JM , 

TEMP) , 

TJM, 

DUMA, 

D’JMC ) 

CALL 

VPLUS 

(  TJF* 

FJET, 

TJF* 

DUMA* 

DUMC) 

955  CONTINUE 

GO  TO  96 
050  CONTINUE 

CALL  S  C  A  L  M  (  TJF,  DV 1 •  TJF,  2.*  DUMA) 
CALL  SCALMITJM,  DV),  TJM,  2.*  DUMA! 
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060 


continue 


r 


r  COMPUTE  MOMENTS  DUE  TO  SECTION  CMS 

r 

DO  961  1=1*  25 

DELM  =DYNPR( 1 )*CMR< I )*CSSW*.04*(CHDK( I >**2) 

DO  961  K  =  1  *  3 

961  TSMIK)  =  TSM(K)  +  DELMMBONDR  (K) 

I F (SYM)  GO  TO  965 
DO  962  1=1*  25 

DELM=  DYNPL(I)  *CML(I)  *CSSW  *.0a  *< CHDL ( I ) ** 2 ) 
DO  962  K  =  1 »  3 

oft?  T  SM ( K I  =  TSM ( K  )  +  DFLM*  ABONDL(K) 

GO  TO  97 

°  6  S  TSM  (2)  =  2.*  TSM<2) 

TSM(l)  =0, 

T  SM ( 3 )  =0. 

970  CONTINUE 

r 

r  FIND  FORCE  AND  MOMENT  COEFFICIENTS 

C 

C  1,  RESOLUTION  TO  STABILITY  AXES 

C 

SAXM)  =-U(  1 ) /SORT  (  U  I  1  )**2  +  U  (  3  )  *  *  2  ) 

S AX ( 2 )  =0, 

SAX(S)  =  -U( 3)/SQRT(U(1  )**2  +  U ( 3 ) ** 2) 

S  A Y ( 1 ) =0 • 

S  A  Y  (  2  )  =  1  • 


SAY ( s ) =0. 
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CALL  VCROSS (SAX, SAY, SAZ, DUMA, DUMB) 

DO  971  K=l,3 

TOTF(K)  =  TF (K)  +  TJF(K) 

971  TOTM(K)  =  TM(K)  +  TJM(K)  +  TSM(K) 
ROXORG(l)  =  XCBAR 
ROTORG(2)  =  0 
ROTORG(3)  =  CBAR 

CALL  VCROSS (ROTORG , TOTF , TEMP , DUMA , DUMB) 
CALL  VPLUS (TOTM. TEMP, TOTM, DUMA, DUMB) 


CALL 

VO0T<  SAZ, 

tote  , 

DV1  f 

CLIFT, 

DUMA  ) 

CALL 

VDOT ( SAX, 

TOTF, 

DV1  > 

CDRAG, 

DUMA  ) 

CALL 

VD0T ( SAY, 

TOTF, 

0V1  * 

CSIDE, 

DUMA  ) 

CALL 

VDOT ( SAZ, 

TOTM  , 

D  V 1  * 

CYAW  , 

DUMA  ) 

CALL 

VDOT ( SAX, 

TOTM, 

D  V 1  ♦ 

CROLL , 

DUMA  ) 

CALL 

VDOT ( SAY, 

TOTM, 

DVlf 

CPITCH, 

DUMA  ) 

CPI  TCH  =  CPI  TCH+TOTF  m*XCBAR  -TOTFd  )*ZC8AR 
C 

r  ADJUST  non-dimens ionali ZI Ng  factors  and  signs 

r 

CLIFT=-CLIFT/ESS 
C DR AG=~CDR AG/ESS 
CP  I T  CH=  CP  ITCH/ (ESS*  CBAR) 

C  S I D  E  =  CSIDE/ESS 
CROLL  =  CROLL /  <  2  *  *  ESS) 

CYAW  =  CY AW / ( 2 • *  ESS  ) 
r 

C  COMPUTE  STATION  LOADINGS  IN  OUTPUT  FORM 

r 

DO  1201  1=1,  2S 

DO  1202  K= 1 »  3 

1202  TEMP(K)=EOR( I »K)  +  FJR(I»K) 

CALL  VDOT(TEMP,  SAX,  DV1,  DLOADR ( I ) ,  DUMC) 

CALL  VDOT(TEMP,  SAZ,  DV1,  'JPLOdR(I),  DUMC) 
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SLOAQR ( I )=TEMP(2> 


UPLODR ( I ) =~UPLODR (  I  ) 

DLOAOR (  I  )=-DLOADR (  I  ) 

5  F  C  C  L  R  (  1  )=  UPLODRI  I  >*25./CH0R(  I  1 

1201  S  FCCDR ( I  ) =  OLOADR t  1  )*25./CHDR<I ) 

IF(SYM)  GO  TO  1210 

00  1203  1=1*  25 

DO  1204  K  =  1  *  3 

1204  TFMP(K)  =FOL ( I ♦  K)  +  FJL(I,K) 

CALL  VOOKTFMP,  SAZ,  DV1*  UPLODL  I  I  )  ,  DUMC) 

CALL  VOOT ( TEMP  *  SAX*  0V1*  OLOAOL(I)*  DUMC) 

OLOADL ( I  )  =  TEMPI 2  ) 

OLOADL ( I  ) =-DLOADL (  I  ) 

UPLODL! I  ) =-UPLODL (  I  ) 

FFCCLL(I)  =  UPLODL  (  I)  *25./a-nL  (  I  ) 

1202  SECCDLII)  =  OLOADL ( 11*25. /CHDL(I) 

1210  CONTINUF 

r 

r  PRINT  RESULTS 

C 

NPAGE=NP AGE+1 

WRITE(6*1001)  DATE1*DATE2»TITLE*  NPAGE 
1 0 0 1  F0RMAT(1H1*  Aio.  A2*27X*10A7/40X»50<1H. ) »22X»*  ^AGE  ** 

112/*  WING  CHARACTERISTICS-  ASPECT  RATIO  SWEEP  ANGLE 

2  TAPER  RATIO  DIHEDRAL  ANGLE  TWIST*) 

WRITF(6*  1002)  aR*  SWFEP*  TaPpR*  OIHdRL.  TWIST 
mo?  FORMAT(]H  *28X,F5 .2 ♦ 1 2X,F5.?**  DEG* » 11 X , F5 . 3 » 10X *F5 . 2  *  *  DEG*. 
U0X,F5,2»*  DFG*/) 
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WRITE(6»1003) < I »  I  =1  » NPAN ) 

1003  FORMAT!*  FLAP  ARRANGEMENT-  *»6<*  PANEL  *»  11)/) 

WR I TE ( 6  »  1 004 )  (  Fp  (  I  )  ,  I  =2  »  NPaNPd 
inn4  FORMAT!*  END  SPAN* , 6X » 6 ( 1 0XF4 . 2 > ) 

WRITE(6» 1005) (OF! I ) *  1=1*  NPAN) 
pnc  FORMAT!*  CHORD  RATIO  *,  6(10XF4.2>  ) 

IF(SYM)  GO  TO  11 00 
W  R  T  T  F ( 6  »  1 006  > ( DFR (  I  )  ♦  1  =  1*  NPAN) 

ion6  pormak*  Right  deflfction  *,6(F6*2**  deg  *)) 

W  R I T  F ( 6 ♦  1007 > ( FFR (  I ) ♦  1  =  1*  NPAN) 

1007  FORMAT!*  WING"  EXTENSION  **  6(F4.2*10X)  ) 

WRITE ( 6, 1008 ) ( DFL (  I  )  ,  I  =1  ,  NPaN  ) 

1008  FORMAT!*  LEFT  DEFLECTION  *,  6<F6.2»*  DEG  *>) 

WRITE (6*  1007 ) ( FFL (  I )  ,1  =  1 ,  NPaN) 

GO  TO  1101 

1100  W  R I T  E ( 6 ♦  1000 > ( DFR ( I ) ,  I  =  1,NPaN) 

1000  FORMAT!*  DEFLECTION*,  12X,  6(F6»2**  DEG  *)) 

W  R I T  E ( 6 »  1 0  1 0  )  (FFR! I )  ,1  =  1*  NPAN) 

1010  FORMAT!*  FXTFNSION  *,  14X,  6(F4«2»  loX)/) 

wRTTF(6,lnll)  (BLOS(I).  1=1,  NPaN) 

1011  FORMAT!*  BLOWING  DISTRIBUTION-*,  7X*  6<F5.3*  OX ) / ) 

GO  TO  1102 

1P1  CONTINUE 

WRITE<6»  1012)  (BLOR(I)*  1  =  1*  NPAN) 

1012  FORMAT!*  BLOWING  D  I  STR  I  BUT  I  ON-*  /  1  5  (  1H  )»*"IGHT  WIiNG  *»6<F5.3*9X 
1  )  > 

WRITE(6»1013)  ( BLOL ( I ) *  1=1,  NPAN) 

101^>  FORMAT  I  1  5  (  1H  ),*LEFT  WING  *»  6  <  F  5  -  3  ,  9X  )  ) 
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no?  CONTINUE 


I r ! SYM )  GO  TO  1030 

WRITE  (6*]031)  (CJPR(D*  1  =  1*  NPAN ) 

ioi]  format ( *  panel  cJ  5-  right  wing  **6< F6.3»8X) ) 

WRlTP ( 6* 1032 )  (rJPLtl),  1=1,  NPAN) 
lOi?  FORMAT ! *  LEFT  WING  *  *  6 < F6 • 3 ♦ 8* >  ) 

GO  TO  1034 
10-^0  CONTINUE 

WRITE(6*1033) ICJPRt I) »  1=1*  NPAN) 
lO'JO  FORMAT!*  PANEL  CO  S-*  *  16X  *  6  (  F6. 3  »  8X  )  ) 
n34  CONTINUE 

I F ( BLOWDR.GE.o. )  WR  I  TE ( 6  *  1 0 14 ) 

IF (BL0WDR.LT. 0. )  WR ITE(6»1015) 

i o i a  format(*  blowing  parallel  to  plane  ^f  symmetry.*) 

1 0 1 5  FORMAT!*  BLOWING  NORMAL  TO  HINGE  LINE#*) 

WR!  Tf<6,i016)  ALff(NaS)*  BETa(NBS),  C-MNJS) 

1016  FORMAT!/*  FLIGHT  CONDITION-  ANGLE  OF  ATTAC*-  SIDESLIP  A^G 

1LE  OVERALL  CJ*/26X»  F6.2**  DEG  *.10X,F6.2»*  DEG* *  1 2X » F5 . 2 ) 

WRITE(6.8001)  NTRY,  EXCESS 
3001  FORMAT!*  TRIAL  *,I3»*»  EX CE SS=* * F 1 0 . 5 ) 

IF (SYM)  GO  TO  1103 

WRITE  !  6  *  1017)  CLIFT,  CDRAG.CSIDE,  CPITCH,  CR'-'LL,  CYAW 

1 0 1 7  FORMAT!/*  FORCE  AND  MOMENT  LIFT  DRAG  SIDE  FORCE 

1  PITCH  ROLL  YAW*/*  COEFFICIENTS** 12X* 

2f6,3,4X,F8.4,  f 10 . 4 , 14XF8. 4 , 5X , F9 , 5  * 3X , F9. 5/ lH  *125(lH.)//*  RIGHT 
OWING  DATA-*/) 

GO  TO  1104 

1103  W  R  I  T  F  !  6  *  1018)  CLIFT,  CDRAG  ,  CPITCH 
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1018  FORMaT(/*  FORCF  and  MOMENT  COEFFICIENTS-  LIFT-*F7.3**  DRAG"*. 
ifp.4.  *  pitch-**  f9.4/ih  »125<ih.i//*  st'ation-by-statiun  data- 

?.*/ ) 

1104  WRITEI6, 1019) 

1.010  FORMAT  { *  STA  SPAN  CHORD  ClRCU-  LOCAL  WIND-  LOCA 

1L  LOADING-  CHORD  SECTION  COEFFICIENTS-*/ 

2*  NO  LOC  (NOM)  LATION  SPEED  ALPHA  BETA  LIFT  ‘  DKA 

3G  SIDE  (EXT)  CL  CD  C J* ) 

DO  1105  1=1.  2E 

wRTTC-(6*1020)  I,  OR  (  I  ,  2  )  *  CNOM(I),  GR  (  I  )  ,  VELR(I).  ALFA*(I).  BETA* 
1(1).  UPLODR ( I  ) .  DLOADR(I).  SLOADR (  I  ) , CHDR ( 1 )  . SECCLR (  I )  . SECCOK ( I  ) 

2,  C JR ( 1 ) 

IF(SYM)  60  TO  1301 

]  0  2  0  F0RMAT(1H  *  I  2 » F6 . 2 ♦ 2F9. 4  * F 1 0 . 3  . F 8 . 2 »  F7.2.  F8.4.2F9.5,  F7.4, 

1 F 1 2 • 3 .  F 1 2 • 5 ♦  F12.3) 

NPAGE  =  NPAGE+1 

WRITE(6.1021)DATE1  ♦  DATE2»  NPAGE 

1021  FORMAT ( 1H1 »  98XA1D.  A2»*  PAGE  *,I2/*  LEFT  WING  DATA-*/ ) 

WRI TE (6. 1019) 

DO  1106  1=1.  25 

1106  WRITF(6. 1020)  I,  0L(I,2)*  CNuM(I),  GL  (  I  )  *  VELL(I),  ALFAM  I  )  .  BETAO 
1(11.  UPLODL ( I  I »  DLOADL(l).  5L0ADL (  I  ) , CHDL (I) , SECCLU  I )  . SECCDL ( I  ) 

2  .CJL(I) 

1301  CONTINUE 

r 

r  SET  UP  NEXT  CASE 


N AS=  NAS+1 
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IF(NAS.LE.NA)  GO  TO  500 


N  A  S  =  1 

NRS=  NBS+1 

I F (NBS.LE.NB )  GO  TO  500 
NRS=1 
N JS=NJ$+1 

IF(NJS.LE.NJ)  GO  TO  A 50 
MJS  =  1 

N EX  T  =  X  I N ( 10) 

GO  TO  (1,  2  *  3*  4*  5 *  5  5*  6*  7*  8.  9,  10*  11*  12.  13)*  NEXT 

CALL  EXIT 

END 

SUBROUTINE  VICTOR  (A,  3*  C.  D*  E) 

DIMENSION  A ( 3 1  *  B(3).  C  £  3  ) 
c.NTRY  VPLUS 
DO  1  J= 1 ♦ 3 
I  C(J)=  A(J)+B(J) 

RETURN 

CNTRY  VCROSS 

Cl  1 )  =  A(2)*B(3)-A(4)*B(2) 

C(?)  =  A ( 3 ) *B ( 1 ) —A ( 1 )#B( 3) 

C  (  3 )  =  A ( 1 ) *B( 2 )“A ( 2 ) *B(  1  ) 

RETURN 
ENTRY  VDOT 

D=4(1)*B(1!+A(2)*B<2)+A(3)#B(3) 

RETURN 
ENTRY  VM AG 

c=A(l)**2  + A ( 2 ) **?  +A(3>**2 

94 


0=  SORT  C  E ) 


RETURN 

ENTRY  SCALM 
00  2  J  =  1  .3 

2  C( J)=  D*  A ( J ) 

RETURN 
ENTRY  VMINE 
00  3  J=  1  » 3 

3  C  (  J  )  =  A  (  J  )  -  B  t  J  ) 

RETURN 

ENTRY  SCALD 
DO  4  J*  I  •  3 
u  S  (  J  )  =  A  (  J  )  /D 

RETURN 
END 

SUBROUTINE  WASH  (P*A*  0,  W) 

LOGICAL  SEGM 

DIMENSION  P(3>*  A<31*  Q  <  3 )  *  W<3)*  R  <  3 »  *  AXR(3)»  Vd<3»*  VDDI3) 
1  *  RR ( 3 ) ♦ rLL ( 3 )  »WD ( 3 ) 

SEGMr. FALSE. 

C0ST=-1. 

in  P l =3, 141 3R3 

CALL  VMINE(G»  P*  R*  X.  Y) 

CALL  VCROSSIA.  R «  A  X  R  ♦  X*  Y) 

CALL  VMAG(  AXR,  VD«  VDD*  X,  nEN ) 

I F (DEN.LT . 1 .OE-lO ) GO  TO  20 
CALL  VMAG (  R*  VD  *  VDD*  RM ,  X) 

I F ! RM.LT.l .0E-10 I  GO  TO  20 
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CALL  VDOT ( A ,  R,  V  D  ,  ADR,  X) 

FAC=( ADR/RM-COST) / (4.*PI*DEN) 

CALL  SCALM ( AXR  »  VD*  W»  FAC*  X) 

I F ( SFGM ) GO  TO  30 

RETURN 

FNTRY  SFG 

^EGM=.TR'JE. 

00  3  K=  T  *  3 
Wn(K)=A(K) 

CALL  VMINE  (0,  A,  RR,  X.  Y) 

CALL  VMINE  (A,  P*  ELL.  X.  Y) 
CALL  VMAGtFLL,  VDD.VD,  AM »  Y) 

I F( AM.GT.l.OE-10)  GO  TO  1 
CONTINUE 
DO  2  K=  1 » 3 
W  (  K  )  =  0  . 

I F ( SE GM ) GO  TO  30 
RETURN 
CONTINUE 
AM=1 ./AM 

CALL  SCALM(FLL,  VD*  A.  AM,  x> 

CALL  VDOTtA,  RR,  VD*  ADRR  *  X) 

CALL  VM AG ( RR  *  VD,  VDD  *  RRM  *  X) 

CO  ST=  ADRR/RRM 

GO  TO  10 

CONTINUE 

DO  31  <=1,3 


A ( <  )  =WD( K ) 


RETURN 


END 

SUBROUTINE  SMOOTH (G) 

COMMON  F ( 2 5 )  ,  T ( 2 5  ) 

0  I  MENS  I  ON  G ( 2  5  >  *  Bill) 

PI =3.141593 
DO  2  K=  1*  11.2 
A  =  0  » 

DO  3  N  =  1 »  25 

8  A  =  A+  F  (  N  )  *G  (  N  )  *  COS(K*T(N)l 

?  B ( K  )  =  4.*A/PI 

DO  4  N=1  »  25 

G  (  N  )  =  0  « 

00  5  K=  1 ,  1  1  .  2 

A  G ( N ) =  G  (  N ) +  B(K>*CO^(K*T(N>  > 

4  CONTINUF 

RETURN 
END 

SUBROUTINE  CEEL ( ALFA .  CJ.  DELTA,  CF  ,  ALFIN,  CL,  CM) 

CM  =  -l.5708*ALFA  -  (  4 . 62*  SQR  T  (  CF  1  -2.93*  CF)*DELTA 
CL=6 .2832  * ( ALFA  +  (.  32+  1  .  155*CF ) *DELTA ) 

I F ( CJ.LE.O. )  GO  TO  1 

CL  =CL+  (2.76*OELTA+  1  .092*ALFA  )  *C  J**.  68  +  (  DEL  T  A  +  A«-FA)*CJ 
F?=1.25*CJ+1.5*(1.  -EXP(  ~1.204*CJ)) 

CM=CM-  DELTA*F2*FXF(-1.18892*CF !  -  ALFA*F2/5. 

1  C  M  =  C  M  +CL/4. 

C  L~  CL  -CJ*SI N( ALF IN ) 


RETURN 
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To  exit,  any  number  ^14.  is  used  for  the  "next  case  card  number". 
If  a  number  less  than  1.  is  used,  an  error  return  will  cause  termination  of 
execution. 


5.  Sample  Problem 

The  externally  blown  flap  problem  discussed  in  Section  IV  is  used 
as  a  sample  problem  here. 

Figure  22  is  the  input  keypunch  form  for  cards  to  run  a  symmetrical 
five-panel  configuration  at  three  angles  of  attack  and  one  Cj,  then  to 
analyze  the  same  case  with  the  fourth  panel  on  the  left  wing  (representing 
the  area  behind  the  left  outboard  engine)  reduced  to  zero  blowing.  (The  flap 
deflection  is  also  changed  to  reflect  the  fact  that  the  flap  centerline  angle 
did  not  coincide  with  the  jet  deflection.) 

Output  for  the  symmetrical  7.5°  angle  of  attack  case  is  shown  on 
Page  99.  The  printed  data  is  mostly  self-explanatory.  The  blowing  distri¬ 
bution  factors  are  given  in  their  normalized  form.  The  "local  loading"  is 
the  force  on  the  strip  divided  by  the  semispan  squared  and  by  the  freestreara 
dynamic  pressure.  It  Includes  both  "circulation"  and  jet  forces.  The  "section 
coefficients"  are  load  per  unit  span  divided  by  local  extended  chord  and  by 
freestream  dynamic  pressure. 

Output  for  the  same  case  with  engine  out  is  shown  on  Pages  100  and  101 
The  same  data  is  given,  plus  the  lateral-directional  coefficients,  the  flap 
arrangement  and  blowing  data  for  the  other  wing,  and  the  station-by-station 
data  for  the  other  wing. 
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